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SYMPOSIUM ON RESERVOIR DATA 


In February 1957 the Kuwait Branch of the Institute of Petroleum organized a symposium 

of papers on the general theme of the acquisition and on the analysis of reservoir data. The 

papers which were presented at this meeting are reproduced in the following pages, together with 
summaries of certain other papers and notes which were made available to the delegates. 


PRELIMINARY CONSIDERATIONS OF ALBORZ OIL RESERVOIR 
AT QUM* 
By H. FARKHAN ¢ and J. TAVANA ¢ 


A cuRSORY examination of the figures representing 
the gross quantity of residual oil lost with the abandon- 
ment of depleted reservoirs * will emphasize the 
importance of scientifically controlling petroleum 
reservoirs and well systems and carefully analysing 
all the fundamental data which may be available. 
The data which may be required have been grouped 


as follows ?: 


(1) Physical properties of the reservoir. 

(2) Physical character of the reservoir fluids. 

(3) Static characteristics of the rock-fluid 
system. 

(4) Dynamic characteristics of the rock-fluid 
system. 


It should be pointed out, however, that reservoir 
studies can neither be restricted to any one specific 
group of factors nor limited to any time or place. 

For instance, it is known that the structural condi- 
tions of the reservoir, and many geological factors, 
may affect its behaviour.*!° Similarly, it may be 
said that the oil reservoir control commences with the 
beginning of drilling operations by controlling and 
adjusting the essential properties of mud in such a 
way as not to plug off the pore system and hamper 
flow of oil,4 etc. Studies of reservoir characteristics 
may be extended outside the reservoir proper, such 
as those which were carried out in the Asmari moun- 
tain and which threw considerable light on the be- 
haviour of Masjid-i-Suleiman reservoir.® 

It was such facts as these which led the Iran Oil 
Company to initiate a series of studies covering all 
the possible characteristics of the prospective Alborz 
oil pool. Attention was first directed to the lime- 
stone formation which is thought to contain commer- 
cial accumulation of oil.!* The rocks considered are 
the limestones which belong to the Marine Formation 
of the Oligo-Miocene * age. 

This limestone was cored in a test well located NE. 


of Alborz and just S. of Separ-i-Rostam (Fig 1), 
where it crops out as an inselberg.4 It was also 
systematically sampled in Kuh-i-Nardaghi (Fig 1), 8. 
of Qum. 

The above-mentioned test well reached the top part 
of the Marine Formation at a depth of nearly 100 m, 
and the formation was continuously cored to a depth 
of 900 m. Cores were recovered in good condition. 
Bulk density, grain density, and porosity measure- 
ments were carried out on samples cut off at almost 
one metre intervals. A summary of these measure- 
ments is tabulated in Table I. 


I 


Porosity and Density Values of Marine Limestone Cores 
South of Separi-i-Rostam 


T 
Grain density | 
range 


Depth range, m Porosity range, 


100-150 
150-200 
200-250 
250-300 
300-350 
350-400 
400-450 
450-500 
500-550 
550-600 
600-650 
650-700 
700-750 
750-800 
800-850 
850-900 


2-53-2-74 
2-48-2-70 
2-51-2-79 
2-56-2-76 
2-53-2-76 
2-55-2-71 
2-53-2-87 
2-58-2-76 
2:60-2-75 
2-56-2-96 
2-59-2-71 
Not available 
2-53-2-73 
2-57-3-04 
2-59-2-73 
2-63-3-01 


These density and porosity values have been 
determined gravimetrically... They were checked 
with values obtained by means of an air expansion 
porosimeter,® as well as by means of the saturation 
method, using kerosine as the saturating medium.™ 

A consideration of the porosity values of the lime- 
stone samples collected at Kuh-i-Nardaghi shows that 


* MS received 17 June 1957. 


+ Technical Director, Iran Oil Co. 


t Head of Petroleum Engineering, Iran Oil Co, 
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6-11-3 
| 1-17 
2-3-20 
4-3-13-5 
3-2-18-4 
| 1-7-16 
| 3-1-18-7 
| 5-7-17-3 
7-24-12-8 
6-47-15-7 
| | 2-04-7-2 
1-73-7-9 
1-37-9-9 
| 2-03—8-7 
1-12-68 
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they vary laterally, increasing with the decrease of 
the carbonate content of the limestone and its gradual 
change into marl. 

The only reservoir fluids available for examination 
were the oils which flowed out of the Alborz No. 3 


Queternory Lekes [EES Juower Red Formetion 
Mio-Pliocene incl, URed. Form. Fj Eocene mainly tuffeceous 
[ZY Selt Ptugs Diepiric zones Igneous Miocene-Eocene 
Form Sligo Miocene) 
Anticlines —*— Synclines 


Strike lines 
soo ke 


Fie 1 
GEOLOGICAL MAP OF QUM AREA 


well and Alborz No. 5 well. The results of the tests 
are tabulated in Tables II and III; the change of 
specific gravity and viscosity of A.5 weathered crude 
are shown graphically on Figs 2 and 3. 

Distillation, viscosity, and asphaltene content tests 
were carried out in accordance with the IP and the 
rest according to ASTM standards. 

As has been mentioned elsewhere,” it was not 
possible to obtain a complete set of data on the wells 
drilled at the Alborz field. The following facts may 
be of interest. 

The upper red formation,* ® overlying the cap rock 
(evaporite) is relatively easily penetrated by the drill. 
However, as soon as the chemical series acting as the 
cap rock is reached and ordinary mud is replaced by 
saturated brine mud to drill through this series, which 
includes rock salt layers of varying thickness, con- 


TaBLE II 
Analysis of Crude Oil from Alborz No. 3 Well at Qum 


Specific gravity at 60° F ; 0-8354 
Redwood I viscosity at 100° F : . 47 sec 


Water and sediment 0-2% vol 
Carbon residue 2-72% wt 
Sulphur 0-13% wt 
Asphaltenes ‘ 0-365% wt 
Chloride content . 11 p.p.m. 
Distillation : 
Distillation to: 
75°C . 
100°C . 
125°C . 
150°C . 10% 
175°C . 16% 
200°C . 21% 
225°C . 26-5% 
250°C . 32% 
278°C . 39% 
300°C . 45% 
Volume of: 
Distillate 46% 
Residue 53-8% 
Loss 0:2% 


TaBLeE III 
Analysis of Crude Oil from Alborz No. 5 Well at Qum 


C2. 

C3. 

Ist-C4 

Normal C4 

Ist-C5 

Normal C5 

C6 and higher . 
Specific gravity 

API gravity . 
1.B.P. on distillation 
Volume of distillate . 


Total distillate 

Total residue . 

Sp. gr. of distillate . 

Sp. gr. of residue 

Total sulphur . 

Carbon residue 
Asphaltenes (hard asphalt) 
Sediment (by extraction) . 


Wax content . 

M.p. of wax 

Reid v.p. 
Viscosity, Redwood I 
Kinematic viscosity 
Pour point (upper) . 
Total acidity . . 
Water . 
B.S. & W. 

Salt content 


Crude oil (drum sample taken at atmospheric pressure) : 


Trace Trace 
0-14 vol % 0-06 wt % 
0-52 vol % 0-32 wt % 
0-26 vol % 0-17 wt % 
1-22 vol % 0:86 wt % 
0-92 vol % 0-69 wt % 
1-39 vol % 1-06 wt % 


95-55 vol % 
0-8310 at 60° F 
38-8° at 60° F 
52°C 

3% at 75°C 
8% at 100°C 
13% at 125°C 
18% at 150°C 
28-5% at 200° C 
33% at 225° C 
38-5% at 250° C 
44% at 275°C 
49% at 300° C 
49-5% by vol 
49-5% by vol 
0-7765 at 60° F (50-65° APT) 
0-8863 at 60° F (28-15° APT) 

0-085% by wt 

1-25% by wt 

0-08% by wt 

0-03% by wt 

0:005% by wt 

il 


Ni 
10-8% by wt 
132° F 


5-8 p.s.i.g. at 100° F 
37 sec at 100° F 

4-57 cS at 100° F 
70°F 

0-024 mg KOH /g 

0-06% by vol 

0-2% by vol 
12-3 Ib/1000 bri 


JOURNAL OF THE INSTITUTE OF PETROLEUM 


| 
% > 


OF ALBORZ OIL RESERVOIR AT QUM 47 


siderable stuck-pipe trouble was encountered. As an 
instance the Alborz No. 4 well may be mentioned. 
The drilling of this well began with a 184-inch hole 
using ordinary mud about 75 lb/cu. ft. and having a 
1500-1000 Marsh funnel viscosity of 32-37 sec. The 
casing, 13%-inches o.d., 68 lb/ft. J.55, was set at 


TEMPERATURE - °F 


80 
SP GR 
0-82 0-83 0:84 0-85 
Fie 2 
170 
1 
ng 
KIN 
Fie 3 


151 m. Thereafter the drilling of the well was con- 
tinued to the chemical series, i.e. to a depth of over 
2000 m with an uncased hole of almost the same depth 
showing no sign of caving and causing no difficulty. 
However, shortly after the sweet water mud was 
replaced by the brine mud, the drill pipe got stuck 
and pieces of rock from the upper part of the long 
uncased hole began falling on the bit. The hole had 
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to be abandoned because the drilling tool, which was 
frozen at a depth of over 2500 m, could not be pulled 
free. 

The next well, Alborz 5, was drilled to the evaporite. 
But before penetrating the chemical series, where the 
use of heavy brine mud is imperative, a second string 
of casing was inserted in the well to protect the upper 
formation from the disintegrating effect of the brine 
infiltration. The abnormal formation pressures of 
the evaporite series, which were considered mainly 
responsible for freezing the drill pipe in Alborz 4, 
were counteracted in Alborz 5 by using a mud of 
128-135 lb/cu. ft. in drilling through that section. 
This time the difficulty of handling a heavy mud, 
such as the hazards caused by the swabbing effect of 
the bit while pulling out, was the main difficulty. 
The Alborz No. 5 well, however, reached the oil hori- 
zon, but only after a series of protracted operations of 
fishing and side-tracking were painstakingly overcome. 
It may be of interest to note that the side-tracked hole 
lies so near one of the fish that the drill often touched 
its end when it was being lowered to the bottom. 

At the depth of about 2400-2500 m, gas was 
occasionally noticed bubbling on the mud surface. 
The testing of the formation was ruled out on account 
of extremely precarious conditions existing at the 
bottom of the hole. No gas bubble, however, was 
reported immediately before the blow-out. 

Thus the drilling of the 84-inch hole went ahead 
with 44-inch drill pipe, 2-5}-inch x 2} id. and 
10-5-inch o.d. x 2-inch i.d. drill collars. At 2676 m 
a change of formation was recorded. At about 
03-20 on 26 August, the well began to surge and rapidly 
unload the 129 lb/cu. ft. mud column and it then 
blew wild. Several streams of oil began flowing away 
from the wild well. The largest one of these was 
gauged and was found to be running at about 80,000 
b.d. Taking into consideration the oil sprayed 
around the well, and the oil flowing in smaller streams, 
it was estimated that the initial oil production 
amounted to about 100,000 b.d. or even more. This 
was borne out by the quantity of oil accumulated in 
several rapidly improvised reservoirs, the construc- 
tion of which involved about 5 million cu. ft. of 
earthwork. Some details of these reservoirs are 
already published.*® 

It should be pointed out that right from the begin- 
ning considerable interest was taken by all concerned 
in evaluating the true reservoir pressure. Thus it 
was estimated that at the moment the well blew wild 
the mud column was exerting a pressure of about 
8000 p.s.i. at a total depth of about 8782 ft (2677 m). 
Gas-cutting could not be contemplated because the 
mud properties were being checked at 15-minute 
intervals throughout the shift. No change either in 
the mud properties or in the position of the mud level 
in the mud pits was observed. 

When the capping assembly was installed the ques- 
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tion of observing the accurate wellhead pressure and 
thereby estimating the reservoir pressure presented 
itself.’ To visualize the situation a brief description 
of the capping assembly may be helpful. The 
assembly consisted from bottom to top of: 


(1) An adaptor spool which was to be connected 
to the 7-inch casing spool, to the southern side 
of which a plug valve had been connected. 

(2) A flow spool with 2-4-inch flanged outlets 
each fitted with two plug valves. When assembly 
was installed one outlet faced E. and the other 
faced W. 

(3) A hydraulically-operated Cameron D.V. 


valve. 


The western flowline was fitted with a 3-inch posi- 
tive choke. The eastern flowline was connected 
through a Shaffer variable choke, having a maximum 
opening of 1 inch, to a 2000-ft flowline. 

A 2-inch W.K.M. gate valve was added to the 
aforesaid 2-inch plug valve on the 7-inch casing spool 
in the wellhead assembly facing 8. One joint of 
23-inch tubing was screwed into the 2-inch W.K.M. 
valve. After the above-mentioned 23-inch tubing 
was screwed into the W.K.M. valve, both valves were 
opened, permitting flow through the tubing. The 
4-inch plug valves on the E. flow line were opened, 
permitting flow through the 1l-inch choke. The 
valves on the W. flowline were closed and the wellhead 
pressure increased immediately to 2500 p.s.i. and re- 
mained constant at this value. 

An attempt was made to close the W.K.M. valve, 
but after two turns of the valve wheel the valve 
could not be closed further, nor could it be reopened. 
The wellhead pressure had, meanwhile, jumped to 4500 
p-s.i. The 4-inch valves were immediately reopened 
on the W. side and the pressure declined to 2500 p.s.i. 
A few minutes later the 23-inch tubing blew out of 
the W.K.M. valve. 

As it was necessary to remove a 4-inch x 2-inch 
reducer, which had previously been installed on the 
end of the 4-inch flowline to the W., the 4-inch valves 
on that side were closed and the pressure built up to 
3500 p.s.i. in the few seconds necessary to make the 
change, and was still increasing when the valves were 
reopened. 

A second attempt to close the 2-inch W.K.M. valve 
was unsuccessful and the 2-inch valve, fitted before it 
on the 7-inch casing spool, was necessarily closed. 
At this time the 1l-inch choke on the E. side plugged 
solidly with debris and cavings from the hole. The 
well was then flowing to the W. through the 3-inch 
i.d. positive choke and 4-inch line. 


In considering the Alborz reservoir pressure it may 
be apposite to point out that the salt plug, W. of the 
over 50-km-long Alborz structure,® is believed to be 
flowing out at a rate sufficient to compensate for the 
continuous removal of the exposed salt by the seasonal 
rainwater. The seat of the high pressure required, 
not only to counterbalance the weight of the Kuh-i- 
Namak salt plug, but also to induce flow at the said 
rate, may be the same as that causing the high pressure 
phenomena observed in the A.5 well. 

These facts, taken in conjunction with the strati- 
graphy * and structural geology of the area, leave 
little doubt as to the extent and the extraordinary 
condition of the Alborz oil reservoir partially mani- 
fested through the A.5 blow-out. The successful 
capping of the A.5 well was due not only to the 
selfless devotion of the IOC rank and file, ably led by 
F. Naficy, Exploration and Production Director of the 
National Iranian Oil Co., and by B. Mostofi, Managing 
Director of the Iran Oil Co., but also to an inter- 
national co-operation, coupled with discipline as well 
as good fellowship, which has often been declared 
impossible by short-sighted sectarianism. 
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WELL INTERFERENCE TESTS IN THE SOUTH DOME OF THE 
ZUBAIR RESERVOIR * 


By M. J. E. HEATON f 


SUMMARY 


Well interference tests were carried out on the three wells in the South Dome of the Zubair Reservoir which 
have been completed in the Third Pay. For this purpose one of the wells (Well No. Zb 17) was put on stream at 
5000 b.d., and minimum closed-in pressure measurements taken at Well No. Zb 2 (3 km away) and Well No. Zb 19 
(10 km away) to determine the form of the resultant pressure drawdowns there. 

Two studies were carried out. One, after 150 days, concentrated mainly on Well No. Zb 2; and the other 


after 200 days, when Well No. Zb 19 was reviewed. 


The theory of point sinks in an infinite reservoir was applied, using the method of images to allow for the 
size of the reservoir in the study of Well No. Zb 19. Selected theoretical drawdown curves were obtained. By 
fitting various such curves to the observed drawdown at Well No. Zb 2 (corrected for overall pressure decline) 
figures were obtained for the reservoir parameters. These figures may then be fitted into the more complicated 


equation for Well No. Zb 19. 


A check is thus available on the parameters, which have also been determined by laboratory investigation. 
The main purpose of the tests, however, is to determine the connectivity of the producing sands, i.e. what per- 
centage of the sand at one well continues to a point some distance away. 


PREVIOUS HISTORY OF THE SOUTH 
DOME 


A rough map of the South Dome, showing the three 
relevant wells, is appended (Fig 1). It will be seen 
that the South Dome is an oval structure roughly 20 
km long by 10 km broad, with the three wells approxi- 
mately in a line down the long axis. The average 


WATER Lever 


NORTH 


912 3 6 5K 
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ROUGH MAP OF SOUTH DOME, SHOWING THE THREE WELLS 
COMPLETED IN THE THIRD PAY 


sand thickness above oil—water level is about 120 
feet. It is not known whether the structure is 
faulted off from the North Dome, but certainly there 
is little if any oil connexion between the two domes. 
The pay is about 50 per cent sand. 

It should be emphasized that these figures are only 
known very approximately, largely on account of the 
lack of data from drilled wells. 

The total production from the South Dome prior 
to the commencement of the tests in mid-May was 


less than 100,000 brl, and there was no production at 
all for the preceding five*months. 

A pressure survey was carried out to determine the 
pressure decline in the dome consequent on the North 
Dome production. Declines of 0-7 p.s.i/month at 
Well No. Zb 19, 1-1 p.s.i/month at Well No. Zb 2, and 
1-7 p.s.i/month at Well No. Zb 17 were found. 
Allowance is made for this effect. 


GEOLOGICAL EVIDENCE 


The Zubair Reservoir consists of a large number of 
thin sands, separated by silts and shales. Owing to 
the wide spacing of the wells, particularly in the 
South Dome, it is impossible to tell with any con- 
fidence whether sands at corresponding depths in 
neighbouring wells are connected or not. A certain 
amount of negative evidence is available, however, in 
that certain sands in some wells obviously have no 
counterpart in their neighbours, or correspond to much 
thinner sands there. 


MATHEMATICAL THEORY 


The differential equation for one well in an infinite 
reservoir is * 
 10P 


where P is pressure drawdown in atmosphere at dist- 
ance r cm from the well at time ¢ sec; 
¢ is porosity, as a fraction; 
k is permeability, in darcies; 
u is oil viscosity, in eP; 
c is compressibility (of oil and sand combined), 
in vol/vol/atmosphere. 


* MS received 17 June 1957. 
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A solution of equation (1), for a well producing at 
constant rate q cc/sec (reservoir oil), from sand of 
thickness h cm, is? 


where Ei is the exponential integral, defined by 

—Ei(--k) =| 
Solutions of the form (2) may be combined 

additively. This allows the method of images to be 
used conveniently to deal with reservoir boundaries, 
faults, etc. 


DISCUSSION 


First we must transform equation (2) to a mote 
convenient system of units: 


4:26 x 10®%ducr? 
Eh kt ). 


where P is pressure drawdown, in p.s.i.; 

r is distance from producing well, in km; 

t is time of production, in days; 

q is production rate, in b.d. (of stock tank oil, 
assuming a formation volume factor of 
1-52); 

k is permeability, in mD; 

h is formation thickness, in ft; 

¢ is porosity, as a percentage; 

u is viscosity, in cP; 

c is compressibility, in vol/vol/p.s.i. 


(4) 


Note that h is average connected sand thickness 
between the two wells in question. 

Fig 2 shows the observed drawdown at Well No. 
Zb 2 after 150 days, corrected for overall pressure 
decline, compared with the curve in 


—aini(— 7). 


which is seen to provide a good fit. 
Remembering that r = 3 for Well No. Zb 2, we have 


3 = 4-26 x 106 x F x de; 


=—1073x% x 
21 = x Xx (6) 

Now laboratory investigations give fairly accurate 
measurements of f and c as 19 per cent and 16 x 10-* 
vol/vol/p.s.i. respectively. Also g = 5000. Hence 
we have 

= 0:00232; h= 593. . . (7) 
Using the laboratory figure » = 0-6, this gives k = 
259, which agrees, within the limits of error, with 
laboratory measurements. 


Turning attention now to well No. Zb 19, and 
assuming that the parameters deduced above remain 
constant (with the exception of the connected sand 
thickness h) we may adopt the general equation 


p= — (-*) 


as the pressure drawdown due to a given well. 
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IMAGE REPRESENTATION OF SOUTH DOME 


The assumption that the parameters other than h 
remain constant will be justified if a good fit is ob- 
tained. As will be seen, any variation in them must 
be small compared to the variation in h. 

From Fig 1, it is clear that one must take into 
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account the size of the reservoir when dealing with 
wells as far apart as Wells No. Zb 17 and Zb 19. 
This is done by the method of images; the field is 
represented by a series of wells in an infinite reservoir, 
as in Fig 3. 


DRAWDOWN PS! 
> 


In practice, it is found sufficient just to use the one 
“image,” at 13 km from Well No. Zb 19, since for 
each of the others the resultant drawdown after 200 
days is of the order of 0-01 p.s.i. or less. This gives 


1250 ,,. 300 1250 ,,. 500 


for the pressure drawdown of Well No. Zb 19 con- 
sequent on the production at Well No. Zb 17. 


for overall pressure decline. Again a good fit is 
obtained. 


CONCLUSION 


The theoretical investigation of well interference 
tests herein carried out gives results in good agree- 
ment with laboratory investigations, and also sheds 
light on the connectivity of the sands in the South 
Dome. 

It is suggested that about half the sand thickness in 
the individual wells connects between Wells No. Zb 
17 and Zb 2, a distance of about 3 km; and about a 
quarter between Wells No. Zb 17 and Zb 19, a distance 
of about 10 km. 
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DISTRIBUTION OF FLUID PROPERTIES IN A HIGH RELIEF OIL- 
FIELD AND ITS EFFECT UPON THE MATERIAL BALANCE * 


By H. K. GRANT 


SUMMARY 


Vertical variation in the fluid properties of high relief oil reservoirs has been ized for several years, but 
there is little or no reference in the Sieustane on its significance in estimating the tank oil originally in place by 
the material balance method. 

Using a hypothetical reservoir saturated at or near the crest with the variation in the fluid properties with 
elevation based upon published data, it was found that no single pv can be used in calculating the oil in place 
by material balances. It was found that the reservoir must be numerically integrated by dividing it into hori- 
zontal layers, each having its own p-v-t characteristics. 

Curves are presented to show the probable variation with elevation of the fluid properties of a hypothetical 
reservoir. Also presented are material balance calculations based upon three methods: (1) using a variable 
p-v—t (numerical integration method); (2) using a volumetrically weighted p-v-t; and (3) using a p-v-t based upon 


the average tank oil produced. 


INTRODUCTION 


THE fact that the physical properties of fluids vary 
with elevation in a single oil reservoir has been recog- 
nized for many years. Sage and Lacey in 1939 did 
some theoretical work predicting changes in the 
physical properties with structural position. Since 
then several authors have noted increasing fluid 
densities with depth, particularly on high relief pools. 

The purpose of this paper is two-fold: first, to 
show, using a hypothetical reservoir, how the fluid 


properties might vary with elevation based upon 
published data; secondly, to calculate the effect these 
variations would have in a material balance type of 
reserve estimate using a volumetrically weighted 
p-v-t and one based upon the average produced oil 
as compared to a more realistic consideration of the 
variation in the p-v—t with elevation. 

It can safely be said that almost all of the work 
done to date on small reservoirs involved the use of a 
single p-v-t. For most of the cases the errors in 
using a single or average p-v-t in estimating oil and 


* MS received 17 June 1957. 
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gas reserves by material balances are probably small, 
but for some the errors could be quite substantial and 
may explain, at least in part, why many material 
balance calculations have given erroneous results. 

The hypothetical reservoir used in showing the 
variation in the physical properties of the fluids and 
the importance of considering the reservoir in hori- 
zontal layers relative to p-v-t data is one having a 
total relief of 1000 ft. Several curves are presented 
showing the variation in the physical properties with 
elevation and pressure. All p—v-t data were deter- 
mined by applying the physical properties to Stan- 
ding’s empirical p—v-t correlations for mixtures of 
California oils and gases.2__ Also presented are example 
calculations for three material balances based upon: 
(1) a variable p-v-t; (2) a volumetrically weighted 
p-v-t; and (3) a p-v-t based upon the average pro- 
duced oil. 


DISTRIBUTION OF FLUID PROPERTIES 


The amount of published data on the estimation 
of oil and gas reserves by material balances and the 
prediction of the future performance of oil reservoirs 
by considering the variation in the physical properties 
of the contained fluids with structural position is 
extremely limited. Probably the first complete reser- 
voir analysis which considered a changing p-v-t with 
elevation is the one presented by D. R. McCord * on 
the performance predictions of the LL-370 area, 
Bolivar Coastal field, Venezuela. It is shown in his 
paper that the oil gravity varied between 18° and 
28° API over 1500 ft of closure for an average gravity 
change of about 0-7° API per 100 ft of elevation. 
In the fifth zone of the Newhall-Potrero field in 
California there occurred a variation of tank oil gravity 
from 25° API to 39° API within a closure of 2000 ft.4 
The average change in API gravity over the range of 
closure in this field was 0-7° API per 100 ft, which 
agrees with the variation found in the Bolivar Coastal 
field. The initial saturation pressure in the Weber 
sand of the Rangely field in Colorado varied from 
2560 to 1960 p.s.i. over 840 ft of relief ° for an average 
saturation pressure change of about 71 p.s.i. per 100 ft 
of elevation. Similar variations have been found in 
the large Burgan field in Kuwait. 

In developing a hypothetical reservoir for the pur- 
pose of illustrating the variation of the physical 
properties and in showing the effect of these variations 
on material balance calculations, an average change 
of 0-7°API per 100 ft and an average change of 70 
p-s.i. per 100 ft of elevation in the initial saturation 
pressure was used to agree in general with known 
data on high relief pools. The variation in the tank 
oil gravity and initial saturation pressure with eleva- 
tion should not be taken as linear. An oilfield with 
a closure‘of 1000 ft between the elevations of 6000 
and 7000 ft is assumed. A circular homogeneous 


reservoir with a 5-mile radius containing a constant 
net sand thickness of 200 ft, a porosity of 28 per cent, 
and water saturation of 10 per cent of the pore volume 
is assumed. Other values for the porosity and water 
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saturation could have been used without affecting 
the conclusions at the end of the paper. 

The change in the tank oil gravity is assumed to 
vary between 28° API and 35° API (Fig 1); the 
reservoir temperature to vary between 140° and 160° F 
(Fig 2). The initial solution ratio was assumed to 
range from 400 to 700 cu. ft. per barrel of tank oil 
(Fig 3). The gravity of the produced gas was assumed 
to vary between 0-77 and 0-83 (Fig 4), and the gravity 
of the free gas was taken as 0-70 over the reservoir 
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conditions used in the material balance calculations. 
Fig 5 shows the original saturation pressure—elevation 
relationship for the assumed reservoir, obtained by 


ery 


6200 


ELEVATION — FEET SUBSEA 


27000! 
700 600 500 400 
SOLUTION GAS-OIL RATIO—CU.FT/BBL 
Fie 3 
GAS SOLUBILITY v. ELEVATION 


applying the physical properties of the fluids in 
Standing’s empirical p-v-t correlations for mixtures 
of California oils and gases. The empirical p—v-t 
relationships established by Standing were obtained 
by running numerous p-v-t on various types of Cali- 
fornia crudes. The relationships are normally pre- 
sented in the form of three charts whereby the 
saturation pressure and liquid phase formation volume 
factor can be obtained by knowing such field data as 
the producing gas-oil ratio of a bubble point liquid, 
and the produced oil and gas gravities plus the reser- 
voir temperature. The variation in the original 


pressure checks the 70 p.s.i. per 100 ft of elevation 
observed in actual reservoirs. In developing the 
original reservoir pressure—elevation relationships also 
shown on Fig 5, it was assumed that the reservoir 
was just saturated at the top. 

The original liquid phase formation volume factor 
was obtained by applying the liquid phase compres- 
sibility obtained from Fig 6 showing the com- 
pressibility above the saturation pressure as a func- 
tion of the saturation pressure. This curve was 
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obtained by drawing a smooth curve through some 
thirty measurements of saturation pressure and liquid 
phase compressibility. It is recognized that this 
curve is only roughly correct and also that the liquid 
phase compressibility above the saturation pressure 
is not constant but decreases slightly with increasing 
pressure. However, this error can be considered 
negligible. Table I shows the calculation of the 
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I 
Calculation of Original Liquid Phase Formation Volume Factor 


| 
Formation 
Reservoir | Saturation Compres- volume 
Elevation,| Gradient, | pressure, pressure, sibility, factor at 
ft subsea | _p.s.i/ft P,, C, per Car Ba; CAPBy 
p.s.i.a. p.s.i.a. 10° p.s.i. Res. Brl/ | 
O8 

Brl T.O. | CAPB,, 
6000 0-3081 2620 2620 13-80 1-360 
6100 0-3095 2651 2605 46 13-73 0-00063 1-354 00009 
6200 0-3105 2682 2590 92 13-65 0-:00126 1-350 0-0017 | 1-848 
6300 0-3120 2713 2565 148 13-53 0-00200 1-343 0:0027 1-340 
6400 0-3140 2744 2535 209 13-39 0-00280 1-335 00037 «1-331 
6500 0-3160 2776 2500 276 13-22 0-00365 1-327 0:0048 | = 1-322 
6600 0-3181 2808 2455 353 13-00 0-00459 1-318 0-0060 1-312 
6700 0-3213 2840 2385 455 12-66 0:00576 1-305 0-0075 1-298 
6800 0-3256 2872 2275 597 12-15 0:00725 1-288 00093 | 1-279 
6900 0-3323 2905 2110 795 11-32 0-00900 1-262 00114 | 1-251 
7000 0-3412 2939 1900 1039 10-30 0-01070 1-230 0-0132 | 1-217 

II 


Calculation of Original Reservoir Pressure Gradient 


= 1 2 3 4 5 6 a 8 9 
Original Initial 
: Oil es Initial | Wt of tank | Wt of pro- — ss Wt of oil liquid reservoir 
Elevation, if it solution | oil per bri, | duced gas oll oil plus gas phase pressure 
ft subsea ne “ patie 0 ratio, 350-16 x | per cu. ft., 0-076362 x per bri T.O.,| formation gradient, 
po lb 0-076362 x XR 44+6, volume | 7 x 0-4331 
Po Po ft/brl Po» Ib Ib factor, |B, x 35016’ 
Bo, bri/brl p.s.i/ft 
§000 0-8498 0-770 700 297-57 0-0588 41-14 338-73 1-360 0-3081 
6100 0-8509 0-771 690 297-95 0-0589 40-64 338-59 1-353 0-3095 
6200 0-8519 0-772 680 298-30 0-0589 40-05 338-35 1-348 0-3105 
6300 0-8534 0-773 665 298-83 0-0590 39-24 338-07 1-340 0-3120 
6400 0-8550 0-775 650 299-39 0-0592 38-48 337-87 1-331 0-3140 
6500 0-8576 0-778 630 300-30 0-0594 37-42 337-72 1-322 0-3160 
6600 0-8602 0-782 607 301-21 0-0597 36-24 337-45 1-312 0-3181 
6700 0-8639 0-787 576 302-50 0-0601 34-62 337-12 1-298 0-3213 
6800 0-8686 0-795 535 304-15 0-0607 32-47 336-62 1-279 0-3256 
2 6900 0-8762 0-808 475 306-81 0-0617 29-31 336-12 1-251 0-3323 
7000 0-8871 0-830 400 310-63 0-0634 25-36 335-99 1-218 0-3412 
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original liquid phase volume factor. The original 
reservoir pressure gradient was checked by calcula- 
ting the density of the liquid and gas phases. This 
calculation is shown on Table IT. 

The relative gas density for the free gas was obtained 
by assuming a free gas gravity of 0-70 at reservoir 
conditions over the pressure range used in the material 
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(5) Reservoir initially saturated at the crest. ae 

(6) All oil is produced between the elevations eS 
of 6400° and 6800 ft subsea. Although this eo 
assumption is valid for material balance case 
No. 1, slight changes in the location of the 
average p—v-t in the upper and lower zones would 
occur, due to drainage into the middle zones. 


PRODUCING YEARS 


Fie 8 


POOL PERFORMANCE OF A GIVEN HYPOTHETICAL RESERVOIR 


balance calculations and determining the com- 
pressibility of the free gas from the gas gravity using 
the relationship that the pseudo critical pressures and 
temperatures are functions of the gas gravity. This 
method is explained by G. G. Brown.* Fig 7 shows 
the relative gas density as a function of pressure. 
Figs 10 and 11 are the working curves for the gas 
solubility and liquid phase formation volume factor 
obtained from the given data using Standing’s p—v-t 
correlations. 


ASSUMPTIONS 


The past performance of the hypothetical reservoir 
and the material balance calculations are based on 
the following assumptions: 


(1) The reservoir consists of four layers, within 
each of which the properties of the fluid are the 
same. This applies only to the first material 
balance case. The layers and mean depths used 
are-— 


6000-6400 ft subsea with mean depth 6300 ft 
subsea. 

6400-6600 ft subsea with mean depth 6500 ft 
subsea. 

6600-6800 ft subsea with mean depth 6700 ft 
subsea. 

6800-7000 ft subsea with mean depth 6900 ft 
subsea. 


(2) Little or no gravity drainage. 

(3) Constant free gas gravity over the range 
of pressures considered. 

(4) Tank oil and produced gas gravities do not 
change with time at a given elevation. 
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However, for this particular case this can be 
ignored due to the relatively low rates’ of oil 
withdrawal. 
(7) Expansion of interstitial water is negligible. 
(8) No form of pressure maintenance. 
(9) No effect from rock compressibility. 
(10) No water drive. 
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MATERIAL BALANCE BY LAYERS 


A standard form of the material balance may be 
written as follows :? 


NAyBo — R) G. 
(Ri — R) — — B,) 


where N = oil originally in place in bri; 
N, = cumulative oil produced in brl = =Nen; 
y = relative gas density « 5-615 cu. ft/brl in 
8.C.F. per reservoir barrel; 
B, = liquid phase formation volume factor in 
reservoir barrels per brl of tank oil; 
G, = cumulative gas production in cu. ft.; 
B,; = initial liquid phase volume factor; 
R = solution gas-oil ratio in cu. ft/brl of tank 
oil; 
R; = initial solution gas-oil ratio. 


N= 


When more than one layer is used in a reservoir to 
consider the variation in the fluid properties the 
material balance equation can be written: 


— R)] + 
fal (Ri — R) y( Bu 


where f denotes the ratio of the tank oil originally in 
place in each layer to the total tank oil originally in 
place in the reservoir. This equation requires the 
introduction of the relative value of the tank oil 
in place from the volumetric estimate and a know- 
ledge of oil production from each layer. 

In showing the comparison of the results obtained 
between p-v-t data that varies with elevation, a 
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volumetric average p—v-t and one based on the average 
produced oil, it is assumed that the tank oil originally 
in place was 12-3 x 10° brl to agree with the volu- 
metric estimate. It should be emphasized, however, 
that other values for the tank oil originally in place 
could be used without affecting the conclusions. In 
other words, the volumetric estimate was made only 
for the purpose of obtaining the f factor used in the 
material balance equation. 

Fig 8 shows the assumed pool performance curve 
and Fig 9 the assumed average pressure—depth rela- 
tionship at production periods of 2, 6, and 10 years. 
The average rate of production is assumed at about 
100,000 b.d. The pressure-depth curves used were 
arbitrarily made to give values of tank oil in place 
of 12-3 x 10° brl, using the material balance case 
where the p—v-t data varies with elevation. In other 
words, recognizing that this type of material balance 
is the most realistic. The other material balances 
are calculated to show how inconsistent their results 
are and how greatly their answers vary from the 
assumed realistic answer. 

Table Ill shows a material balance calculation 
based upon a producing period of ten years, using the 
method of dividing the reservoir into four separate 
layers. From the calculations it is clearly evident the 
most critical factor in the material balance is the 
saturation pressure-elevation relationship which 
determines the value of f(R; — R). Because of the 
critical nature of this factor, the reservoir is separated 
into small layers, each 100 ft in thickness, in the 
region where the pressures have fallen below the 
initial saturation pressures. 


III 
Material Balance on a Given Hypothetical Reservoir after a Producing Period of Ten Years using a Changing = 
| | ‘Top to | 6400-6600 | 6600-6800 | 6800-7000 | 
— R)] + |6400fts.s.,| ft s.s., ft s.s., ft s.s., 
N = R) — B | Symbol | Unit 6300 ft | 6500 ft | 6700 ft | 6900 ft Total 
| | pvt, p—v-t, p-v-t, p-v-t, 
| f = 0-1098 | f = 0-1915 | f = 0-2917| f = 0-4070 
Cumulative production: | 
Gas = <4 M.c.f. x | 129 129 
Pressure. ‘ P.s.i.a. 2078 2105 2172 2322 
Relative gas density ‘ 8.C.F/R.C.F. 157 158 163 166 
phase fv B, _| Res. brl/brl T.O. 1-275 1-272 1-276 1-258 
Solution gas—oil ratio . . and R 8.C.F/brl T.O. — 510 507 — — 
| §.C.F/brl T.O. one 619 661 
NA5615yB, — R)/1000 Muc.f. x 108 118 126 
N5-615yB, — R)/1000 + G, | Mue.f. x 247 | 255 502 
Initial phase f.v. Bx Res. brl/brl T.O 1-340 1322 | 1-298 1-251 
(Bu — | Res. brl/brl T.O.| + 0-065 | + 0-050 | + 0-022 — 0-007 
y(Ba — B.) + 10-205 + 7-900 | + 3-586 — 1-162 
5-516y(Bui — Bo) F 8.C.F/brl T.O + 57-301 + 44-359 | + 20-135 — 6-525 
—f[5-615y(Bo — S.C.F/brl T.O. | — 6-292 | — 8-495 | — 5-873 +2-656 | —18-004 
—R) . 8.C.F/brl T.O. | 16-624 23-416 18-612 58-652 
| | 40-648 
Tank oil in place T.O.LP. Brl x 10° | — 12-34 
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MATERIAL BALANCE USING WEIGHTED 
AVERAGE p-v-t 


The weighted average p-v-t was based upon the 
mean gravity of the tank oil originally in place. The 
mean elevation of this average oil is taken from Fig 12 


+ 
i 
be T T 

4 4 +. + 

t 
= + 
+ 
4 
- OR ORIGINALLY IN PLACE — PERCENT 
Fie 12 


PER CENT TANK OIL ORIGINALLY IN PLACE U. ELEVATION 


For a homogeneous circular reservoir having a 5-mile 
radius with constant dip and uniform net productive sand 
thickness of 200 ft. 


at about 6750 ft subsea. A material balance after 
ten years of production using the p-v-t at this eleva- 
tion is shown in Table IV. In analysing the calcula- 


TaBLe IV 


Material Balance on a Given Hypothetical Reservoir after a 
Producing Period of Ten Years 


A = using a weighted average p-v-t 
B = using an average tank oil p-v-t 


7 A B 
Nm Ge _ | symbol Unit 6750 ft | 6600 ft 
(R, — RB) — y(Boi — Bo) p-w-t p-v-t 
Cumulative production: 
Oil . Ne Brl x 10* 380 380 
Gas Ge M.c 10° 258 258 
Pressure 2200 2132 
Relative gas density 8. ol F/R. 6. F. 165 161 
Bo Res. bri/br1 T.0. 1-279 1-274 
211-035 | 205-110 
Solution gas-oil ratio . ‘ R 8.0.F;brl T.O. | 512 507 
— 8.C.F/br1 T.0. | 673 645 
Be 'R)/1000 M.c.f. x 10* 256 245 
615yB, — R)/1000 + 
° M.c.f. x 10 514 504 
Initial’ liquid phase fv. f Boi | Res. bri/bri T.O. 1-288 1-312 
(Boa — B, Res. bri/briT.0.| +0-009 | +0-038 
y(Bo +1:485 | +6-118 
— — Be) 8.C.F/br1 T.0. | —8-337 | —34-353 
(R, — R). ‘ — 43 123 
Tank oil in place . TTOLP. Brl x 10° 14-83 5-68 


tions it is clearly evident that this type of p-v-t does 
not consider the average effect of the reservoir pressure 
on the solution ratio but only the effect at one eleva- 
tion. Although this material balance gave a value 
of tank oil in place of 14-8 x 10° brl with an error 
of only 16-9 per cent, the material balance after six 
years of production using this p-v-t resulted in a tank 
oil in place of 52-1 x 10° brl with an error of over 
300 per cent. 


MATERIAL BALANCE USING A p-v-t 
AT AN ELEVATION OF 6600 FT SUBSEA 


A material balance based upon the average eleva- 
tion of produced oil after the ten-year production 
period is shown in Table IV. The calculated tank oil 
in place in this case was only 5-7 x 10° brl or 54 per 
cent lower than the given value of 12-3 x 10° bri. 
In this case too much weight was placed on the change 
in the solution ratio at this elevation and time relative 
to the average change in the reservoir. 


SUMMARY AND CONCLUSIONS 


Table V shows a summary of the results of the 
three material balance methods at two, six-, and 


TABLE V 
Summary of Material Balance Results 
Method 
Weighted) Average 
Changing average | tank oil 
p-v-t 
Assumed tank oil originally 
in place, billion brl 12:3 12:3 12:3 
Calculated tank oil in place 
(two years’ 12-3 16-6 18-4 
Error, % + 35:0 | 
Calculated tank oil in 1 place 
(six years’ 12-3 52-1 6-6 
Error, % +3236 | —463 
Calculated tank oil in place 
(ten years’ 12-3 14:8 5-7 
Error, % . ‘ -—~ + 169 | --53-7 


ten-year periods and how they differ from the given 
value of the tank oil in place. 

Separating the reservoir into layers for the purpose 
of a material balance becomes even more important 
when an initial gas cap is present. 

Partial water drive could have been assumed for 
the hypothetical reservoir, but the comparative results 
would have been quite similar to those obtained. In 
a given highly undersaturated reservoir it becomes 
evident that the results of a material balance is depen- 
dent entirely on the compressibility of the liquid 
phase, and therefore the proper p-v-t would be one 
based upon the produced oil as long as the reservoir 
is everywhere above the saturation pressure. 

The following are the important conclusions that 
can be derived from this paper: 


(1) The most realistic method of determining 
the value of the tank oil in place by the material 
balance in an oil reservoir that is saturated at or 
near the crest is one that considers the variation 
in the fluid properties with elevation. This is 
clearly evident from the inconsistency of the 
other two material balance results. 
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(2) Where the fluid properties vary with eleva- 
tion in a reservoir saturated at or near the crest, 
the saturation pressure-elevation relationship is 
the most critical factor in estimating the reserves 
by the material balance method. 
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APPENDIX 


The following development is to show that there is very 
little error in the material balance results due to oil migration 
by assuming that all of the production is obtained from the 
interval between the elevations 6400 ft subsea and 6800 ft 
subsea 
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The gas balance for the assumed production block may be 
written as: 


[NB; — (N-N.)B — Na Buly + (1) 
or: 


NR; + — ¥Bm) = 
(W-NJR + (NB, — 


N = oil originally in place; 
R; = original solution gas—oil ratio; 
N, = cumulative tank oil; 
R = solution ae ratio; 
Nm = migrated oil 
R,, = solution gas-oil ratio for oil; 
B; = initial liquid phase formation volume factor; 
B = formation volume factor; 
Bm = formation volume factor of migrated oil; 


(N-N.)Bly + @ (2) 


where: 


y = relative gas density for the free gas in the 
reservoir; 
G, = cumulative gas produced. 
The maximum error in assuming that the migrated oil term 


Nun(Rn — yBm) is equal to N,,(R — yB) can be obtained by 
comparing the maximum and minimum probable values for 
N,,, R,, and B,,. 

It can be assumed that in calculating the oil reserves by the 
material balance the maximum probable oil that will be 
drained into the producing block would not exceed 25 per cent 
of the tank oil originally in place in the reservoir when the 
cumulative produced tank oil is about 50 per cent of the oil 
originally in place. 

Further the maximum probable difference between 

NnlRn — yBm) and N,,(R — yB) for the reservoir conditions 
used in the report is obtained by assuming that the migrated 
oil is all from below. This condition, however, can never be 
attained because there will always be some oil migration into 
the producing zone from above. The maximum possible 
value for (Rm — yBm) — (R — yB) will not exceed 22 per cent 
of Ry. 
Therefore the maximum error introduced by neglecting the 
second term on the left side of equation (2) is 0°25 x 22 = 5°5 
percent. This is also the error in the calculated value of N, 
the tank oil in place. 

In general the material balance will only be me Say when 
the cumulative production N, is very much less than 50 per 
cent of N and N,, is proportionally small. Moreover, part of 
the drainage will be up and part down so that the difference 
in the values of yB — R will be less. Hence the error intro- 
duced by neglecting the term for the migrated oil will be far 
less than 5:5 per cent. 
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ANALYSIS OF PRESSURE DECLINE * 
By H. SAMARA + 


SUMMARY 


A well-known principle, the principle of superposition of states, is used as a basis for a new approach to 


the analysis of pressure decline curves. 
production. 


This principle is used explicitly or implicitly by every analysis of variable 
In the present investigation no other principle is used explicitly, and the analysis is valid for any 


laws of flow and for any geometrical boundaries of the reservoir, as long as the principle of superposition remains 


true. 


The analysis enables the prediction of: (1) the bottom hole pressure differential from short flow in wells; (2) 
future pressure decline in the reservoir, provided that the reservoir has been produced long enough. 


EXPERIMENTAL PROCEDURE 


(1) Bottom hole pressure differentials are usually 
measured in the followingway. The well is left closed- 
in long enough for a static bottom hole pressure to be 
established, the pressure is measured, the well is 
allowed to flow at a constant rate until the bottom 
hole pressure stabilizes, then the pressure is measured, 
the well closed in, and, after the pressure stabilizes, 
the bottom hole pressure is measured again to ensure 
its return to its original pressure. The difference 
between the static and flowing pressures is the bottom 
hole differential pressure for that rate of flow. 

The only drawbacks to this procedure are: 


(a) the time of flow for stabilizing bottom hole 
pressure may be unduly long, and this is especi- 
ally inconvenient for a new well; 

(b) since the flowing bottom hole equilibrium 
pressure is not known, the only criterion for 
judging whether equilibrium has been reached is 
the rate of change of pressure with time. In 
certain cases this may be too small to notice long 
before equilibrium pressure is reached. 


(2) In the producing life of a field, pressure decline 
is measured periodically in observation wells, and a 
record is kept of production. 

A plot of pressure drop v. time shows not one smooth 
monotonic curve which could be easily analysed, but 
a curve with ups and downs due to changes of produc- 
tion rates. A pressure decline v. production rate 
reflects changes, but to a lesser degree. These 
changes make the prediction of the pressure behaviour 
of the reservoir by graphical extrapolation rather 
uncertain. The usual procedure is to postulate a 
physical model of the reservoir, work out the behaviour 
of the reservoir for constant production, and deduce 
from that the behaviour for the prescribed production 
and, finally, to match the pressure history with the 
observed one. When a good match is obtained, an 


extrapolation is carried out into the future based on 
the physical model. 

The reliability of the prediction depends to a 
large extent on the length of the record of observation. 

The method proposed here would enable one to 
produce a pure pressure decline curve, i.e. a curve for 
constant production rate, from an actually observed 
curve with varying production rates. If the reservoir 
has been produced for a period long enough to obtain 
the full pure pressure decline curve, i.e. for the 
initial transient pressure drop to have died away, 
then this curve can be used directly for the purposes 
of pressure decline prediction, thus eliminating the 
necessity of introducing a physical model. Even 
when such a model is required, it is easier to deduce it 
from the pure curve than from the observed curves. 


THEORETICAL ANALYSIS 


Principle of Superposition of States 

The principle states that the effect of production on 
pressure decline is additive, ¢.e. that the pressure 
decline at any given time due to a production rate 
which is the sum of several production rates is the 
sum for the same given time of the individual pressure 
declines due to each production rate separately. 
The mathematical expression of the principle is as 
follows. 

If, when a reservoir is produced at any rate q,(t) its 
pressure decline at time ¢ is F(t), and when it is 
produced at rate q,(t) its pressure decline is F,(t), 
then when it is produced at production rate q,(t) + 
q,(t) its pressure decline would be F(t) + F,(t). 

(1) This principle depends on the linearity of the 
equations of motion, which is the case for the Darcy 
single fluid flow equations. However, if the equa- 
tions of motions are non-linear, they can usually be 
replaced for a first approximation by linear equations, 
and the smaller the rates of flow under consideration 
the better the approximation. 


* MS received 17 June 1957. 


+ Mathematician, Petroleum Engineering Dept, Iraq Petroleum Co. Ltd. 
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(2) It also depends on the fact that the conditions 
of the reservoir, like permeability, porosity, and its 
physical boundary conditions, are unaffected by the 
flow. 

(3) The principle does not depend on the geometrical 
boundaries of the reservoir, or on the distribution of 
porosity or permeability inside it. 

(4) The principle implies that the B.H.P.D. is 
proportional to production rate, which is not true in 
fractured limestone fields iri the Middle East. 

Nevertheless, the principle is implicitly used when- 
ever a model of the reservoir is employed to deduce 
pressure drop for variable production rates. 

(5) Let the pressure decline for unit constant pro- 
duction rate be f(t), where ¢ is the time since produc- 
tion commenced, then the pressure decline for con- 
stant production rate q is gf(t). If production started 
at ¢ = 7’, and remained constant at the rate g, then 
the pressure decline is zero from t=0 to t=T 
(since there has been no production in this period) 
and is gf(t— 7) fort >T7. Or since f(t) = 0 for 
t < 0, the pressure decline for this case is gf(t — T'). 

The function 


(I) 


for 7’, more than 7’, represents the pressure decline 
due to a constant production rate g from time 7’, to 
T,. This is so because the first term represents the 
pressure decline due to production rate q beginning at 
time 7',, and the second term represents the effect of 
production —q beginning at time 7',; and hence the 
combination represents the effect of production g from 
T,toT,. It is, of course, zero for ¢ less than 7’. 

(6) The above result can be generalized for the step 
function of production 


q, from 0 to t, 

dn , to te 
to e the pressure decline 


= — — + — 4) — f(t — 
(f(t — — f(t —t,)]+ - 


(7) For the case of a continuous production rate 
q(t), the above equation becomes 


D= f(t—ujdu. . (3) 


(8) The pure curve, D = qf(t), which is obtained 
for a constant production rate g, can be further 


analysed into the following 
D = afit) = qlat + 
. . (4) 
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where a and 6} are constants and F(t) tends to 6 and 
G(t) tends to zero as ¢ tends to infinity. To obtain 
the physical significance of this substitution, put the 
values (4) in (1) with 7, = 0; 7, = T to obtain 


D=4qaT + 6) 


when ¢ tends to infinity, D tends to gaT’ = Qa, where 
Q is the total production. Hence a is the pressure 
decline per unit reservoir production. 

From equation (4), when ¢ becomes large, D tends 
to q(at +6). Since at is the permanent pressure 
drop in the reservoir, 6 is the proper bottom hole 
pressure differential for unit production rate. 

G(t) is thus the transient caused by the start of the 
production. 


Bottom Hole Pressure Differential 

When the bottom hole pressure differential is 
measured the flow is usually short enough to make the 
permanent pressure decline agt due to production 
negligible. Hence from equation (4) 


D=qb—G@t)) ... ©) 


The bottom hole pressuré differential gb is the limit. 
of D when t becomes large enough to make G(t) 
negligible. However, this time may in certain cases 
be very long, and in any case it may be difficult to 
judge when G(t) in fact became negligible, since the 
entity measured is D not G(t). 

However, if in (5) one finds D,, Dr.,, ete., and 
at the same time neglect a, one obtains for 0 << t<T 


D, = — 


Da-ypr+t = ovr + t)— + t) 
Hence 


D+ Dir +... + 
= qG(0) — qG((n — 1)T +2) 


(7) 
But from (6), qG(0) = qb = bottom hole pressure 
differential. 

Hence the bottom hole pressure differential is 
obtainable from equation (7). A special case of 
equation (7) is when ¢ = 7’, when (7) becomes 

Bottom hole pressure differential 

If the term at is not neglected, one obtains 
Bottom hole pressure differential = g@(0) 

= (D,—D) + —D) + — D) + 


where D is limit of D as t tends to infinity. 
In this particular application this equation is only 
of academic interest. 


In the limit 
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The experimental interpretation of measuring 
bottom hole pressure differential using equation (7} 
or (7’) is the following: 


(1) Measure static bottom hole pressure when 
well is closed in. 

(2) Flow for a short duration, making sure that 
the flow rate is constant. 

(3) Measure the bottom hole flowing pressure 
at time ¢ after start of flow; and hence find 
pressure drop D,. 

(4) Close in well at time 7' after start of flow, 
(where 7’ is more or equal to t), preferably soon 
after measurement of flowing bottom hole 
pressure. 

(5) Measure closed-in bottom hole pressures at 
times ¢ + 7',t + 27... after start of flow, and 


differential is assumed to be unity and the time scale 
is arbitrary. The pressure build-up curves have been 
calculated for different duration of flow 7’. In fact, 
curves OAa, OBb, OCc, ODd, OEe correspond to 
T = 0-2, 0-4, 0-6, 0-8, 1-0 respectively. The following 
table shows D,, (the pressure drop at time 7’) for 
different values of 7 and also the sum of Ds for each 


| 
T Dy Ds | Ds, Sum Ds 
O-2 0- 77229 16700 | | 004371 | 0-01250 | 0-00334 0-00116 | 1-00000 
0-93929 | 0-05621 |0-0045 | — 1-00000 
0-6 0-983 0-017 | — 1-60000 
O-8 0-9955 00065 | — — ome 1-00000 


The last column illustrates that the sum is in- 
dependent of the duration of flow and equal to the 
bottom hole pressure differential. 
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thus find the pressure drops D,, 7, D;, 27, - 
ete. Keep measuring the pressure drop at equal 
intervals 7' until the bottom hole pressure returns 
to its original value before the start of the 
experiment. 

(6) Sum all the Ds to obtain the bottom hole 
pressure differential for the given rate of flow. 

(7) In the academic case where the pressure 
comes back to a different value, and thus the 
reservoir undergoes a finite pressure drop D during 
the experiment, subtract D from each D measured, 
and add those to obtain the actual bottom hole 
pressure differential. 


As a graphical illustration of the above, the pressure 
build-down for a fractured porous medium has been 
calculated and plotted in Fig 1. This is the curve 
OABCDEF, where the total bottom hole pressure 


Analysis of Pressure Decline Curves: Special Case 


Divide the past production history into equal 
intervals of time 7 such that in each interval the 
production rate is constant. 

Let the production rates be 


q, from 0 to T 
te from T to 2T 
and In rt to nT. 
Then the pressure decline, (see equation (2)) is 


= alfit) — ft — T)) 


+ anlft yr) — — nT)} 
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Remembering that f(t) = 0 for ¢ less than or equal 
to 0, one obtains 


D, = “f(T) 
D, = alf(2T) — f(T)) + f(T) 
Ds = alf(3T) — f(2T)) + gf f(2T) — + asf (T) 


= alf(nT) — fin — + fin — — 
S(n — 2T)) +... + 


when D, is the pressure decline at time r7’. 
Using the abbreviation 


the above equations reduce to 


D,= “fi 

Dy = fi + 

Ds = + + Ufs 

Dy = Ubi + he + + 


If the Ds and gs are known accurately, then f, is 
obtained by solving the first equation, this value is 
substituted in the second equation, and the value of 
f, is found, etc., until finally f, is found. Once all 
the f,s are found, and since f(0) = 0, equation (8) 
gives the values of f(r7’) for r = 1, 2,. . . i.e. the 
pure pressure decline curve is obtained. When the 
reservoir has been produced long enough, the transient 
in f(t) (see equation (2)) would have vanished and f, 
becomes a constant with respect to n. Thus f, 
becomes known for all n, and equation (9) can be 
used to predict the pressure drop for any production 
history directly, provided its production rates can be 
treated as constant during the unit intervals 7’. 
Otherwise, f, can be plotted against , and the values 
of f, be found graphically or otherwise and these 
inserted in (2) to find the appropriate pressure drops. 

A model of the reservoir is thus unnecessary for 
purposes of prediction. However, when such a model 
becomes necessary, the pure curve f(t) can be used to 
determine the constants in the model of the reservoir. 

The difficulty in applying equation (9) lies in the 
inaccuracies inherent in observation of the pressure 
drop D and rate of production g. A method of trial 
and error, which may need the services of a digital 
computer, may be necessary to find a series of regu- 
larly varying f, that tend towards a limit which best 
fit equation (9). 

Another method of approach is to substitute for f, 
the value (cf equation (4)) 


af + F, 
where F, = F(nT) — F((n — 1)T) 
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Thus equations (9) become 

D, =a, + 

D, = aQ, + + 

Ds = + WF; + + 

Dy = + + +. (10) 
where Q, is the total production up to time r7’. 

The part aQ represents the permanent pressure 
decline in the reservoir due to the production Q, and 


the rest is the transient due to variation of production 
rates. 


Analysis of Pressure Decline: General Case 


When the restriction that the production rate is the 
stepwise function referred to in the previous para- 
graph is removed the pressure decline p is 


p=fat—wfwdu . ay 


deducible from equation (3) by the substitution v = 
t—u; and the pure decline curve is the solution of 
this integral equation. Integrating (11) by parts one 
obtains 


P = qoyf(t) + [ q(t —u)f(ujdu . (12) 
Hence 


S(t) q(o) q(o) [a (t u)f(ujdu 


or 


S(t) = g(t) + u)f(u)du . . (18) 


P —q (t — u) 
where g,(t) and k(t,w) are (0) and q(o) respect- 
ively. The solution of equation (13), which is equi- 
valent to (11), is possible. As a matter of fact the 
infinite series 
S(t) = + gilt) + galt)... (14) 
where 


Gn(t) = [ . . (15) 


is convergent and satisfies the integral equation (13), 
and thus is the required solution. To prove the 
convergence of (14) it is noticed that both & and g, 
are bounded, i.e. there exist numbers M, N such that 


| k| <M and | g,| 


hence from (15) 


G2 


4 


ig 
| | | 
= 


64 SAMARA: ANALYSIS OF PRESSURE DECLINE 


and it can be proved by induction from the above 
equation that 


| gn] < NM al 
and thus Ifl<=|gn| < 


but the right hand side is convergent (exponential 
series), hence the series for f,, is absolutely convergent. 

Equation (15) can be transformed by integrating 
by parts of the recurrence relation 


which is a more convenient form for computation. 
When the production rate is constant, k = 0, and 
hence g, = o for m more than zero, and thus 


which is what one expects. 

The drawback to this analysis is that there is an 
over-emphasis on the initial production rate q(0), and 
if that happens to be in error, this error is carried 
through the calculation, since it occurs in the denom- 
inator in both equations (15) and (16). On the other 
hand, one can look on q(0) as a parameter, and obtain 
solution (14) for different values of this parameter and 
choose that solution which gives the best fit when 
substituted back in equation (11). A thorough 
investigation of this point deserves notice and shall 
be attempted in the near future. An important 
point in favour of this solution is that one knows that 
it is convergent. 

More general solutions of the integral equation (11) 
can be obtained which are free from the above draw- 
back of the over-emphasis on the initial production 
rate, but which are conditional; i.e. they are solutions 
only in so far that they are convergent. Consider 
any weighing function Z(t), and the series of function 
G,,(t), where 


and G,(t) = Z(t)p. 
Then if the series 


fit) (18) 


is convergent, and thus the function f(t) exists, then 
this function is a solution of the integral equation 


(11). To prove that, add the equation (17) for n = 1, 
2,. .., to obtain 


t 


or 
=f— [at — 

But G, = Zp 

hence G, = Zp = af ae — u)f'(u)du 


which is equation (11) multiplied by Z. 

So far, Z is any function of ¢ which ensures the 
convergence of (19), but from general considerations 
its choice should be such as to make the convergence 
as fast as possible, and thus when the production rate 
is constant q, G, should reduce to p/q, or the value of 


Z(t) for constant production should be ; 


Among the possible functions Z that satisfy this 
general condition the two simplest are 


Z= dZ= 
Fo and = 
The first choice of Z gives us the solution (14) which 
was discussed above. The second choice makes (17) 
read 


G,(t) = G,_,(t) — — u)@’,_,(u)du. (20) 


In effect here i instead of being the initial produc- 


tion rate becomes the average production rate since 
the start of production. 


CONCLUSION AND SUGGESTIONS 


The analysis of pressure decline in general leads to 
two applications: (1) the estimation of bottom hole 
differential pressures; and (2) the construction of a 
pressure decline curve for constant production rate 
from an actually observed curve for varying produc- 
tion rates, and thus to predicting pressure decline 
for any given production. The analysis, though 
theoretical in nature, was general and depended on 
very few plausible assumptions, and thus should have 
wide applicability. 

(1) As far as the bottom hole pressure differential 
is concerned it is suggested that experiments should 
be carried out to test the truth of the theory and to 
establish a routine for such measurements. A series 
of experiments should be carried out in which the 
bottom hole pressure differential, D, of a well is 
measured for a constant rate of flow, q, using the 
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procedure outlined above where the well is allowed to 
flow for time 7. The experiment is repeated for the 
same q and different 7’. The theory shows that the 
bottom hole pressure differential should not depend on 
T, as long asqisthe same. Practically, there may be 
some dependence for small 7’ because of the effects 
of after production, and of the difficulty of keeping q 
constant. The result of these experiments would 
show what is the least practical value of 7 which can 
be used to give consistent results. 

(2) The analysis of pressure decline curves by this 
method presents the arithmetical problem of solving 


equation (9) or (15) and (16). It is suggested that a 
thorough investigation be carried out, possibly with 
the help of a digital computer, into these computa- 
tional problems, especially on the effect of small 
variations of the measured experimental data on the 
final curve. 
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DENSITY OF RESERVOIR CRUDE OIL, WITH REFERENCE TO THE 
CALCULATION OF BOTTOM HOLE PRESSURES * 


By D. E. BUNYAN,? T. I. SHAFIQ,+ and G. J. M. MANASTERSKI f 


SUMMARY 


By measurement of the density of reservoir oil under the conditions found in a closed-in well, bottom hole 
pressures may be calculated accurately from measurements of surface pressures (minimum closed-in pressures). 
The determinations of density are carried out with the aid of conventional p-v-t equipment, with the addition 


of an extra pressure-cell for weighing. 


Weighings are done on a bullion balance. 


The accuracy of the method 


is discussed, and it is shown that the probable error in the pressure computed at the bottom of a 3000-ft well is 


+0-25 p.s.i. or less. 


To ascertain the pressure in an oil reservoir, it may 
often be simpler to measure, in a static well, either 
the pressure at the surface or, in a very low pressure 
field, the level of the free oil, than it is to employ any 
of the various devices for measuring sub-surface 
pressures. The surface measurement is not of value, 
however, unless it can be used for deriving a reason- 
ably accurate sub-surface pressure, and for this pur- 
pose the density of the oil in the well must be known, 
under the conditions prevailing there. W. J. Baker ¢ 
has given a solution to the problem, and the work to 
be described here is largely an adaptation of his work. 
The method of measuring the density itself will be 
discussed, and also the method of employing the 
results in obtaining bottom hole pressures. First, 
however, the method used in determining the surface 
pressure must be considered briefly, in order to evalu- 
ate what range of conditions the density measure- 
ments must embrace. For this, the practices at the 
Kirkuk field will serve as illustration. 

The surface pressure is measured on a closed-in well, 
after the column of oil in the well has been conditioned 
by bleeding off gas as it is formed. A strict routine is 
observed in this procedure, to ensure that the con- 
ditioning will as nearly as possible be the same on 
successive occasions. The surface pressure (or the 


level of the free oil in the open well) is then measured 
at a fixed time (24 hr + 1 hr) after the initial closing- 
in. The result is stated as a gauge pressure in pounds 
per square inch, at a common reference level for all 
wells; in the conversion from the level of measure- 
ment to the level of reference, a pressure-gradient in 
oil of 0-348 p.s.i/ft is assumed, based on earlier but 
somewhat dubious measurements. The pressure so 
computed is called the minimum closed-in pressure 
(M.C.I.P.). More reliance is placed on measurements 
made at a producing well, shut down only for the 
conditioning and the measurement, than on those at 
an observation well, where prolonged shut-down may 
age the oil in ways difficult to predict. 

The pressure, temperature, and degree of saturation 
of the oil in the column therefore range from one set 
of limiting values at the surface to another set deep 
in the reservoir, and the density varies accordingly. 
Kirkuk values will serve as examples. In tempera- 
ture, the range is from a mean value of about 80° F, 
just below the surface of the earth, to something 
between 100° and 120° F in the reservoir. The pres- 
sure will range between zero, at the top of a closed-in, 
bled-off well, to around 1500 p.s.i.g. The saturation 
pressure in the oil near the top of a conditioned well 
will be the same as the ambient pressure, and in the 
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reservoir itself will be between 550 p.s.i.g. and just 
under 400 p.s.i.g. Hence, to be relevant to a study of 
the Kirkuk reservoir, the measurements of density 
must cover a range of conditions from atmospheric 
pressure (at zero saturation pressure and temperature 
of 80° F) to pressures of 1500 p.s.i.g. (at saturation 
pressures between 350 and 550 p.s.i.g. and tempera- 
tures between 100° and 120° F). 

Saturation pressure can be varied by collecting 
unsaturated oil from a point deep in the reservoir and 
bleeding off gas, in the laboratory, until the degree of 
degassing required is achieved. Alternatively, oil 
can be collected from a static, bled-off well, from 
points at various heights above the saturation level 
so selected as to yield oils of the required range of 
saturation pressures. The second course was the one 
followed by Baker, and also in the current set of 
experiments, while the first course was followed in 
earlier measurements made in Kirkuk. The differ- 
ences between the two sets of results will be indicated 
later. 

In his measurements Baker used the sample-taker 
in which the oil had been collected as his “* p—v-t cell.” 
The earlier measurements of density made in Kirkuk 
followed the same method, but the technique is cum- 
bersome. The pressure can be independently varied 
only by injecting water into the sample-taker, and 
the volume changes determined only by weighing. 
The frequent making and breaking of connexions 
makes for many leaks, and weighing is a slow and 
tedious method of finding volume. Standard p—v-t 
equipment and methods were therefore used in later 
work, with certain changes made necessary by the 
special requirements of the experiment. 

Samples were collected in a flush-through type of 
sample-taker (the kind made by the Coleman Instru- 
ment Manufacturing Company) with a capacity of 
750ce. Forthe main part of the experiment, a number 
of samples were collected from various depths in a 
single well on the same day, within a few hours either 
way of the 24th hour from the initial closing-in—the 
well having been given the standard conditioning as 
for M.C.I.P. measurements. The well chosen was a 
producer, the one on which the best continuous series 
of measurements had been made in the past, and to 
which measurements on other wells are correlated. 
The depths of sampling were selected so as to give 
a good range of saturation pressures. Additional 
samples had to be taken from other wells, to obtain 
oil of higher saturation pressure than that now 
obtained from the “ reference ”’ well (but which was 
needed for evaluating earlier readings on the well, 
made when the saturation pressure was higher). No 
sample could be taken from near the surface of the 
oil, as the pressure there is zero and the valves in the 
sample-taker will not close satisfactorily, nor remain 
closed, without some excess of pressure on the inside. 
Instead, a sample was collected from a little below 
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the surface, where the oil was gassy enough to main- 
tain a pressure in the sample-taker, under the temper- 
ature of storage, of around 25 p.s.i.g. This sample 
was subsequently degassed in the laboratory, at 
90° F, to atmospheric pressure. 

The actual measurements were carried out in 
standard p-v-t equipment—pressure cell, constant- 
temperature bath, and mercury pump—manufac- 
tured by Engineering Products Inc. (with pressure- 
gauges by Heise), plus an additional pressure-cell 
made in our own workshops. The samples as col- 
lected were stored in the original sample-takers until 
required. On each, the density was measured at 
1000 p.s.i.g. and 90° F, the saturation pressure was 
measured at three temperatures (one of them being 
90° F), and the volumetric behaviour at the three 
temperatures and over a range of pressures was 
investigated. With the density known at the refer- 
ence conditions of 1000 p.s.i.g/90° F, and the volume 
known at these and a range of other conditions, the 
density at the other conditions is readily computed. 
The procedure in dealing with the successive samples 
changed somewhat as the technique was improved, 
but the latest is more or less as described below. 

From the oil in the sample-taker, some 50-100 cc 
are transferred to the p-—v cell. The saturation 
pressure at one of the three temperatures is measured, 
and this fraction then thrown away. 

The extra pressure-cell, which for brevity will be 
denoted as the ‘“‘ density cell,”’ is thoroughly cleaned, 
dried, evacuated, and weighed on the bullion balance. 
About 250 ce of oil are transferred to the p-w cell, and 
brought to a temperature of 90° F and a pressure of 
1000 p.s.i.g. The density cell is then connected with 
a short length of tubing to the p-v cell, and density 
cell and tubing both evacuated. With the valve on 
the density cell closed, the tubing is filled with oil 
from the p-v cell, which is still kept at the standard 
conditions, and the reading of the mercury pump is 
noted. About 200 ce of oil are now pumped from 
p-v cell to density cell, and the reading of the mercury 
pump again noted. From the difference between 
initial and final readings can be found the volume at 
1000 p.s.i.g. and 90° F of the oil which went into the 
density cell (and is, of course, no longer under those 
conditions), of which the mass will not have been 
affected by the operation. The density cell is then 
cleaned and dried externally, and again weighed; 
the difference between this and the initial weight gives 
the mass of oil in the cell. From this mass, and the 
volume occupied by this oil at the reference conditions, 
the density of the oil at those conditions is computed. 

The oil in the density cell is discarded. Of the oil 
remaining in the p-v cell, the saturation pressure is 
measured at one temperature, and the oil is discarded. 

A further 250 ce of oil are then transferred to the 
p-v cell, and the volume at 1000 p.s.i.g. and 90° F is 
found. The volume is then measured under all the 
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conditions of temperature and pressure that are of 
interest in the particular case—the three temperatures 
at which the saturation pressure has or will be meas- 
ured, and four or five pressures at each temperature. 
The sample is then returned to 1000 p.s.i.g/90° F 
(which provides the opportunity of measuring the 
volume under reference conditions if this was not 
done initially), and the transfer of 200 ce from p—v 
cell to density cell is made, for a check on the value of 
density that was first found. 

The oil in the density cell is discarded. The 
saturation pressure of the oil (50 cc) remaining in the 
p-¥ cell is measured at the third of the three tempera- 
tures of interest, and this oil, too, is discarded. 

If any oil remains in the sample-taker, it is trans- 
ferred to the p-v cell and the saturation pressure is 
checked, at one of the temperatures previously 
measured. 

In all the operations, care is taken not to let the 
pressure of the oil drop below the saturation pressure, 
except in the actual measurements of this property. 
This is because, with this equipment, oil and gas 
cannot be recombined, after having been allowed to 
separate, without large errors in volume being 
observed. The errors are believed to be caused by 
the trapping of gas in the oil—mercury emulsion which 
is formed when the p-v cell containing oil and mercury 
is shaken. 

The observations are repeated on each sample in 
turn, and so the density of the oil under a wide range 
of temperatures and pressures is found, for oils of 
varying gas content. The results of this are given 
in Table I. Before discussing the results, some 
consideration should be given to the sources of error 
in the measurements. The smallness of the changes 
in density produced by small changes in temperature, 
pressure, and saturation pressure, together with the 
closeness with which these functions can be controlled 
or measured (+0-1°F, +2 p.s.i.), makes these 
relatively small sources of error. The volume can be 
measured with an error of around +0-005 cc, which 
again contributes but little to the total error. The 
largest source of error is therefore the weighing. 

The weight of the oil is found as a difference between 
the filled weight and the empty weight of the density 
cell. The cell being necessarily somewhat heavy— 
in this case it weighed just over 6 kg—only a large, 
rather insensitive balance can be used, such as the Oert- 
ling bullion balance which was employed in this work. 
As the oil will weigh something like 200 g or less, great 
care is required in the use of the balance, in order not 
to sacrifice accuracy too seriously. The performance 
of the balance was therefore checked by weighing the 
weights against each other in a cumulative series of 
measurements which found out the sensitivity of the 
balance, assessed the change in balance-point with 
increase in load, determined the repeatability of the 
measurements as given statistically by the standard 
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deviation of a series of readings, and also calibrated 
the weights in terms of one weight chosen as standard. 

A shift in balance-point with increase in load was 
found, and compensated by double-weighings, with 
the cell and the weights switched to opposite pans 
between the two weighings of each pair. 

Calibrating the weights did not sufficiently reduce 
the uncertainty in the values of the larger weights. 
Therefore care was taken to employ the same large 
weights in all weighings. The weight of the oil, 
found as a difference, was then a function only of the 
values of the smaller, more accurately known weights. 

The sensitivity of the balance was found to be 
0-12 g per division of the pointer scale. As the 
equilibrium position of the pointer can be estimated 
by the free-swing method to within one-tenth of a 
division, the error in a single balancing is therefore 
+001 g. Statistical evaluation of the errors in a 
series of weighings gave a standard deviation about 
one-third of the above. As the final weight of the oil 
is the difference between two double weighings, the 
error in this weight is of the order of +0-015 g or 
+0-02g. The contribution of such accidental factors 
as dirt on the cell, or in the fittings on the cell, does 
not seem to have increased the error much above this 
value. 

Taking the weighing error as +0-02 g, and taking 
the weight of the oil as, say, 160 g, then for an oil 
density of around 0-8 g/cc the error in the density 
will be +0-0001 g/cc. Allowing the other sources 
of error to increase this slightly, consider the error in 
density as 0-00012 g/cc. The authors are interested, 
eventually, in computing a bottom hole pressure; if 
the error occurred in the same sense through the 
whole series of measurements (an unlikely event), 
then the calculations of a bottom hole pressure from 
an assumed accurate surface pressure would, for a 
column of 3000 ft, be in error by +0-16 p.s.i. Of 
course, in practice, uncertainties as to the temperature 
and saturation pressure of the oil in the well will 
increase the error, but even so, the method would 
appear to permit the calculation of the pressure at the 
bottom of a 3000-ft well with a probable error of not 
more than +0-25 p.s.i. 

Refer again to Table I, where the density is tabulated 
as a function of saturation pressure, temperature, and 
pressure. For use in calculations of pressure, the 
data must be smoothed, which can be done either 
graphically or algebraically; it was found convenient 
to fit an empirical algebraic equation to the data. A 
reasonably good fit was found for the following 
equation : 


= 0-82816 + 0-03697e— 902552 (4-106 + 
0-000471s) x +- (3-68 + 0-00113s + 0-0147t) 
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where ¢ is the density of the oil in g/cc; 
8 is the saturation pressure (p.s.i.g.) as measured 
at 90° F; 
p is the applied pressure, in p.s.i.g. ; 
t is the temperature, in ° F; 


(The last term in Equation (1) can be ignored at 
saturation pressures greater than 100 p.s.i.g.) 

The last column of Table I shows the deviations of 
the experimental values from the values computed by 
Equation (1). It will be seen that the errors tend to 
be of the same sense and magnitude within the set of 
readings for any given saturation pressure, #.¢. for any 
particular sample. This is what would be expected 
if either the weighing or the saturation-pressure 
determination were the largest source of error; the 
trend of previous argument was that the weighing 
was probably at fault. The standard deviation 
computed for the whole assemblage is +0-000107 
g/cc. If the two values which display errors four and 
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Density at 1000 p.s.i.g. and 90°F as a function of saturation 
pressure 
Graph A = Equation (1) 


five times this are rejected, according to normal 
practice, the standard deviation then works out at 
+0-000074 g/cc. Either value agrees well enough 
with the error anticipated on the basis of the accuracy 
in weighing. 

In Fig 1 the density at 1000 p.s.i.g. applied pressure 
and 90°F is graphed as a function of saturation 
pressure, the experimental values being shown as 
circled dots, and the graph of Equation (1) being 
drawn as a solid line. Also shown are values from an 
earlier, incomplete set of measurements on oils taken 
from various depths in an observation well, from which 
no flow of oil had occurred for several years. The 
third set of values shown are for a sample which was 
taken from the region of the reservoir where satura- 
tion pressure is high, and which was then degassed 
in the laboratory, the density being measured at 
each stage. In these latter two sets the density at 
the highest saturation pressure falls on the line of 
Equation (1), showing that the measurements prior 


TaBLeE I 
Satura- | 
tion Errors, 
pressure | Tem- | Density | g/cc 
at pera- at Pressure, | Volume, | Density, x 105 
ambient | ture, 1000/90 | p.s.i.g. cc g/ee 
ture, 
p.s.i.g. + - 
0 | 80 qa) 526 | 201-833 | 0-83732 | 19 — 
0-8359 356 201-982 | 0-83670 | — 2 
(2) 200 | 202-171 | 0-83592 | — — 
0-8359 52 | 202-470 | 0-83468 | 52 — 
0 90 (3) 1000 202-192 | 0-83583 1 
0-8357 533 | 202-746 | 0-83355 | — 4 
| 359 | 202-957 | 0-83267 | — 3 
198 203-160 | 0-83185 1 
50 | 203-360 | 0-83103 7 — 
| 100 1000 | 203-139 | 0-83193 | — 4 
| 526 | 203-713 | 0-82959 | — 14 
| 357 203-932 | 0-8287 12 
| 201 204-135 | 0-82787 | — 9 
54 204-374 | 0-82691 | 11 
28 80 (1) 1000 | 272-018 | 0-83714 | — 7 
0-8330 464 | 272-859 | 0-83456 | — 11 
} 352 | 273-038 | 0-83402 | — 12 
| (2) 254 | 273-199 | 0-83352 | — 10 
0-8331 153 | 273-234 | 0-83342 | — 49 
32 {| 90 1000 | 273-354 | 0-83305 5 — 
| 464 | 274-223 | 0-83041 | — 1 
352 | 274-406 | 0-82986 | — 3 
251 274-578 | 0-82934 | — 2 
| 148 | 274-751 | 0-82882 | — 1 
35 | 100 1000 | 274-640 | 0-82915 | — 2 
463 | 275-542 | 0-82644 | — 10 
355 | 275-729 | 0-82588 | — 10 
252 | 275-912 | 0-82533 | — 8 
152 | 276-097 | 0-82478 | — 5 
120 80 ql) 1000 | 274-980 | 0-82677 | — 2 
0-8227 499 | 275-813 | 0-82427 | — 2 
402 | 275-983 | 0-82376 | — ‘eon 
(2) 302 | 276-159 | 0-82324 2 — 
0-8228 202 | 276-265 | 0-82292 | — 16 
125 90 (3) 1000 | 276-317 | 0-82277 | — 4 
| 0-8228 530 | 277-125 | 0-82037 | — 6 
410 | 277-348 | 0-81971 | — 2 
299 | 277-551 | 0-81911 1 aon 
200 | 277-733 | 0-81857 4 — 
130 100 1000 | 277-665 | 0-81877 | — 6 
528 | 278-529 | 0-81623 | — 2 
415 | 278-744 | 0-81560 2 — 
303 | 278-955 | 0-81499 3 _ 
204 | 279-138 | 0-81445 5 i 
197 80 999 | 273-583 | 0-82136 6 _ 
605 | 274-251 | 0-81936 5 —_ 
501 274-438 | 0-81880 8 — 
899 | 274-616 | 0-81827 9 -— 
302 | 274-795 | 0-81773 | 14 
207 | 90 1) 999 | 274-943 | 0-81729 7 — 
0-8172 605 | 275-643 | 0-81522 8 ~ 
411 276-000 | 0-81416 | 12 _ 
2) 297 | 276-214 | 0-81353 | 15 —_ 
0-8174 256 | 276-293 | 0-81330 | 17 _ 
218 100 999 | 276-315 | 0-81323 7 eis 
602 | 277-050 | 0-81108 8 — 
500 | 277-239 | 0-81052 | 10 _ 
404 | 277-422 | 0-80999 | 11 ~ 
306 | 277-615 | 0-80943 | 14 _ 
452 90 (1) 996 | 257-396 | 0-80648 | — 4 
0-80662 850 | 257-652 | 0-80568 | — 5 
(2) 698 257-921 | 0-80484 | — 4 
0-80635 595 | 258-106 | 0-80426 | — 3 
478 105 996 | 259-402 | 0-80024 | — 5 
844 | 259-686 | 0-79937 | — 6 
700 | 259-954 | 0-79854 | — 6 
592 | 260-163 | 0-79790 | — 4 
505 120 996 | 261-445 | 0-79399 | — 6 
854 | 261-726 | 0-79313 | — 5 
700 | 262-029 | 0-79222 | — 5 
601 | 262-231 | 0-79161 | — 3 
524 90 (1) 1147 | 268-786 | 0-80504 8 — 
0-8042 1007 | 269-039 | 0-80428 5 _ 
(2) 867 | 269-310 | 0-80347 8 _ 
0-8043 702 | 269-626 | 0-80253 | 10 _ 
553 105 1147. | 270-886 | 0-79880 4 ~ 
| 1006 | 271-163 | 0-79798 4 _ 
| 851 | 271-464 | 0-79710 2 — 
702 | 271-757 | 0-79624 1 _ 
583 120 1148 | 273-026 | 0-79254 3 —~ 
1005 | 273-321 | 0-79166 5 _ 
847 | 273-646 | 0-79074 1 _ 
710 | 273-930 | 0-78992 | — 
| | 
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to degassing are consistent with the present set. 
The density at any particular saturation pressure is 
therefore very dependent on the method of degassing, 
as is shown by the way the sets of values diverge at 
the lower pressures. 

Having obtained an expression for the density of 
oil under a range of conditions which includes those 
in a well conditioned for measurement of minimum 
closed-in pressure, we need only a little further informa- 
tion in order to compute the bottom-hole pressure 
from the M.C.I.P. The temperature in the well 
must be measured, also with the well conditioned as 
for an observation of M.C.1.P., at a sufficient number 


DS/DT - PSI/°F 


100 200 300 400 300 00 
S = SATURATION PRESSURE AT (00°F -PSIG 
Fie 2 


RATE OF CHANGE OF SATURATION PRESSURE WITH TEMPERA- 
TURE AS A FUNCTION OF SATURATION PRESSURE AT 100°F 


DS/DT = 1-95-4-24 (650-8)? x 10-* 


of depths that a reasonably accurate temperature— 
depth curve can be drawn. The saturation-pressure 
of the oil produced by the well must be known, and 
also, of course, the M.C.I.P. 

The variation of saturation-pressure with tempera- 
ture must be studied. The evidence of the present 
study as well as of earlier work is that, to a sufficient 
degree of accuracy, the saturation pressure as measured 
varies linearly with the temperature of measurement, 
but that the slope of this variation is a function of 
saturation pressure. This function takes the form 
shown in Fig 2, where results of a number of investiga- 
tions are assembled. The line drawn fits this equa- 
tion: 

ds 


= 1-95 — 4-24 (650 — s)? x = (2) 


where ? is the derivative of saturation pressure with 


respect to temperature (p.s.i/°F) ; 

sis the saturation pressure (p.s.i.g.) measured 
at 100°F; 

t is the temperature in °F. 


The line is not fitted by statistical methods to the 
whole set of points, but is drawn through what seemed 
on mainly subjective grounds to be the most reliable 
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ones. The scatter being so great, the method is as 
reliable as the data. 

We are now in a position to calculate the bottom 
hole pressure. We take the M.C.I.P. and remove the 
correction to common reference level, to obtain the 
pressure at the top of the column of oil in the particular 
well. From the temperature measurements, pre- 
sumably smoothed algebraically, we have the tempera- 
ture at this point. The saturation pressure at ambient 


« 
Ss \ PRESSURE 
2000+ \ 
=x 
» 
a 
a 
2500F \ 
‘ 
0348 0-350 0352 0-354 
PRESSURE GRADIENT — PSI/FT 
490 690 890 1090 1200 1400 
PRESSURE — PSIG 
Fie 3 
GRAPH OF PRESSURE AND PRESSURE GRADIENT AGAINST 
DEPTH 


M.C.LP. = 155 p.s.i.g. 
Saturation pressure = 476 p.s.i.g. at 100°F 
457°8 p.s.i.g. at 90°F 


temperature is equal to the pressure itself in this, the 
saturated region of the column, and the saturation 
pressure at reference temperature (90° F) is obtained 
from it by use of either Fig 2 or Equation (2). From 
the pressure, temperature, and saturation pressure, 
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the density of the oil is computed, using Equation (1), 
and, from this, the pressure gradient at the top of the 
oil column. With the gradient so found, we compute 
the pressure at a convenient interval of depth (say, 
200 ft) farther down, and by a similar set of calcula- 
tions find the gradient here. The mean of the two 
gradients is then used to find a corrected value for the 
pressure at the lower depth; ifashort enough interval 
is selected it will not be necessary to recompute the 
gradient at the lower point, as the error involved in 
retaining the first-fonnd value is small. With the 
corrected pressure as the new starting-point, the 
whole procedure is repeated, to obtain the pressure 
at a still lower depth. The repetition is continued 
until, with a final depth-interval that will probably 
be of odd length, the saturation level is reached, the 
depth at which the pressure equals the saturation 
pressure (at ambient temperature) of oil from the 
reservoir itself. From here on down, it must be 
assumed that the saturation pressure remains con- 
stant (at least in the Kirkuk case). The density will 
vary—with normal temperature gradients—much 
less with depth; much larger intervals can be used 
in the iterative calculation, chosen to coincide with 
intervals of more or less constant temperature 
gradient. Eventually, the required pressure, that at 
the sub-surface reference level (2200 ft below sea-level 
in the case of the Kirkuk field) is obtained. In the 
course of the calculation, a table like that in Table II 
will have been built up. 

The results of the calculation are shown graphically 
in Fig 3. 

One would not, of course, wish to carry through a 
calculation like the foregoing for every conversion of 
a surface to a bottom hole pressure. The difference 
between the two will in fact vary extremely slowly— 
the change has, according to calculations made as 
above, been only 4 Ib in the whole history of the 
Kirkuk field. It is therefore sufficient to make the 
computation for a few values of the M.C.I.P., with 
appropriate values for the saturation pressure, and 
construct an interpolation table for rapid calculation. 

Calculation of bottom hole pressures by the method 


just described gives results consistent with measure- 
ments by bottom hole gauges. As the calculation 
method should be somewhat more accurate than the 
other, a fairly careful and cunning set of tests would 
have to be made if direct measurements were to 


TABLE IT 
Cor- 
Depth Tem- | Satura- Cor- rected 
below a. Pressure,| pera- tion rected | satura- aoe 
it p.s.i.g. ture, | pressure,) pressure, | pressure, ft 
ft oF p.s.i.g. | p.s.i@. tion 
p.s.i.g. 
0 155-00 86°72 | 157-98 0-35439 
200-00 
200 225-88 | 88-34 | 227-94 225-68 | 227-74 | 0-35239 
200-00 
400 296-16 89-96 | 296-21 295-99 | 296-04 | 0-35069 
200-00 
600 366-13 91:58 | 363-53 365-99 | 363-39 | 0-34929 
200-00 
800 435-85 93°21. | 430-13 435°72 | 430-00 | 0-34802 
83-84 
883-84 464-88 93-89 | 457-80 464-88 | 457-80 | 0-34755 
483-66 
1367-5 632-98 97-81 | 457-80 632-90 | 457-80 | 0-34724 
785-00 
2152-5 905-48 | 101-15 | 457-80 905°50 | 457:80 | 0-34729 
1196-00 | | 
3348-5 | 1320-86 | 105-04 457-80 | 1321-06 | 457-80 | 0-34763 


M.C.LP, = 155-0 p.s.i.g. 

Saturation pressure = 476 p-8.i.g. at 100° F 

Saturation pressure = 457°8 p.s.i.g. at 90° F 
Pressure at ba” ft below sea-level = 1321-06 p.s.i.g. 
B.H.P. — M.C.1.P. = 1166-06 p.s.i.g. 


provide a check on the calculations. Since comple- 
tion of the density measurements, opportunities of 
making tests on producing wells have been somewhat 
limited. It is hoped, however, that the results can 
be corroborated by bottom hole measurements in 
due time. Meanwhile it can be stated that the method 
gives every evidence of being a good one for obtaining 
values for bottom hole pressures from surface 
measurements. 
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SUMMARIES OF TECHNICAL NOTES PRESENTED AT THE 
KUWAIT SYMPOSIUM 


Some OPERATING DiFFICULTIES EXPERIENCED WITH THE 
Ruska SUBSURFACE SAMPLER. P.D.O’Connell. (Qatar 
Petroleum Co. Ltd.) 


The Ruska subsurface sampler is of the flush type. The 
valves of the sample container are “‘ O-ring ” sealed, and are 
closed by the pressure of the well fluid after the timing clock 
has released the tripping mechanism (Fig 1). The sample is 
transferred to the shipping bottle by a closed circuit after 
mercury has been pumped into the shipping bottle to raise the 
pressure to that of the sample-taker, thereby opening the 
valves and allowing displacement of the sample by gravity. 


Difficulties have arisen in the operation of the sampler. 
The commonest trouble has been the obtaining of samples 
with a low opening pressure when connected to the shipping 
bottle, pointing to leakage past a valve. This was overcome 
by machining the ‘‘ O-ring’ groove deeper, and by making 
the diameter of the part of the valve above the ring slightly less 
than that of the part below, so that the lower part of the 
valve seats first. 

Pipe scale and other solids sometimes prevented valves 
from seating properly. A second filter fitted to the upper 
valve gave satisfactory results. Some non-representative 
samples were the consequence of premature operation of the 
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tripping mechanism due to vibration. Careful running of the 
instrument removed this trouble when the sampler was new, 
but weakening of the hammer lock spring on continued use 
made the device more prone to premature tripping. It is 
hoped that springs of Ni-Span C will reduce this weakening, 
which is thought to be a temperature effect. 


OlL CHAMBER S 


CLOCK HEAD 


UPPER VALVE TIMING CLOCK 


CLOCK CHAMBER 
SAMPLE BARREL is 


TRIP MAGNET 


i: Ain LOWER VALVE LOCK \ 
LOWER VALVE S= 
HAMMER LOCK 
HAMMER 
— SCREEN AIR CHAMBER 
RUPTURE PIN 
Fie 1 


In the transfer operation the lower valve did not open fully, 
and the upper valve tended to close under the transfer flow. 
These troubles were overcome by removing the closure spring 
of the lower valve before fitting the transfer head, and by 
installing a stronger spring in the upper transfer head. 

The various modifications improved the reliability of the 
instrument, and led to the obtaining of satisfactory sub- 
surface samples in nine cases out of ten. 
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Dritt-stem Trestrne Resuutts. P. D. O'Connell. (Qatar 
Petroleum Co. Ltd.) 


Conventional Johnston tester equipment has been used in 
Qatar, generally comprising packer, and equilizing, retaining, 
and trip valves, with other auxiliary components as required. 
Single hook-wall packers have usually been employed in 
casing tests, and dual wall packers in open-hole tests. Water 
cushions have ranged from zero to a few thousand feet. 
Open-hole anchors up to 100 ft long have been made up from 
drill collars. Normal packer setting times have been one or 
two hours; however, as much as 30 hours has been used in 
open-hole and 55 hours in casing. Weight on packers has 
varied from 13,000 to 90,000 lb. There has been no trouble 
from stuck packers. Pressure recorders have rarely been 
employed, but the Johnston recorder has given more reliable 
service than the Amerada recorder. 

Casing packer tests have almost always been mechanically 
successful, whereas nearly 50 per cent of the open-hole tests 
have been failures. Most of the failures have been with the 
packers, and have been attributed to excess weight. 


A Nore on THE Use or Bryacat O1m-BasE Mup IN THE 
DvuKHAN Fretp. Qatar Petroleum Co. Ltd. 


The No. 3 and the No. 4limestones of the Dukhan field have 
been cored in a few wells with Binacal oil-base mud in order to 
obtain cores for accurate connate-water determinations. 80- 
brl batches of mud were prepared on the rig site by mixing 
47 bri of gas oil, 8-4 tons of Binfical (a coarse powder consisting 
of a mixture of blown bitumen, calcium naphthenate com- 
pound, naphthenic acids, calcium oxide, and chalk), and 4 bri 
of fresh water. This gave 75-lb/cu. ft. mud. Weighting with 
barytes or calcium carbonate was possible. 

The mud was very stable. When necessary oil was added 
to check rises in weight, and Bitca (blown bitumen and 
calcium oxide) was added to compensate for the diluting 
effects of oil. Contamination with water caused a rise in 
viscosity, but experience and restrictions in the use of water 
on the rig enabled this trouble to be avoided. The initial mix 
had a viscosity of 85-100 sec (Marsh), and an initial gel 
strength of 2-4 g and a 10-min gel strength of 5-8 g; the 
filtrate was zero. 

Drilling with 16,000—-18,000 Ib on an 8§-inch bit, the rate of 
penetration in the anhydrite was about 1-5 ft/hr (half the rate 
obtained with water-base muds). A 5}-inch diamond core 
barrel penetrated the limestone at 2 ft/hr, being slower than 
with water-base mud, but core recovery was 100 per cent in 
nearly all cases. 

The cores were sealed in tins to avoid evaporation losses. 
Crushed samples from the cores were extracted with toluene 
in a modified Dean and Stark apparatus. Consistent results 
were obtained for the connate water content, even though at 
times the water content of the oil-base mud was as high as 20 
per cent. 


ExampLe oF INTERPRETATION OF RADIOACTIVITY AND 
Resistiviry Logs my a SHaty Sanp. P. A. Poupon. 
(Schlumberger International.) 


The earlier empirical formule developed for the interpreta- 
tion of electric logs were applicable only to formations con- 
taining little clayey matter. Later methods were proposed 
which could be used when there was a moderate shale content, 
provided that mud filtrate penetration was at least 4 inches, 
that the residual hydrocarbon content of the invaded zone was 
low, and that the electrofiltration component of the 8S.P. was 
low. Some other information was needed in addition. 

In the Burgan series, however, the invasion is often very 
small and/or the residual oil saturation is high, while S.P. 


NI 
NN 
N\ 
= 
— PISTON af 
Ws, 
4 YA 
NL N 
| 


72 SUMMARIES OF TECHNICAL NOTES PRESENTED AT THE KUWAIT SYMPOSIUM 


deflections opposite clean formations are only 50—60 per cent 
of the value predicted by the conventional formula. As a 
result it does not seem possible to use the Poupon—Loy— 
Tixier interpretation charts. Gamma ray and neutron logs 
have therefore been employed. 

The gamma ray log is basically a recording of the shale 
content of the rocks, provided that the shales have a constant 
amount of radioactive material, and that such material occurs 
only in the shales. Ideally, if a sufficient number of cores 
were available, and their shale content could readily be 
obtained, it would be possible to ascertain whether there was 
a simple relationship between shale content and gamma ray 
deflections. Failing this, it seems that in some cases a 
relationship between shale content and gamma ray deflections 
can be obtained from the logs. 

For the Third Burgan sand, which was water-bearing in a 
Magwa well, a plot of neutron against gamma ray deflections 
gave points which were fairly close to the average curve, and 
it appears that the dominant physical parameter in the rela- 
tionship is the shale content, because of the parallelism 
between gamma ray and S.P. curves. When the average 
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curve was replotted as the neutron-derived porosity index 
against the gamma ray deflections there was an almost linear 
relationship. Correction for the high sodium chloride content 
of the water in the interval studied yielded a relationship 
shown by the broken curve of Fig 2. This suggests the prob- 
ability that the gamma ray deflection and the porosity index 
are approximately linear functions of the percentage of shale. 

The porosity index in a shaly formation is the sum of the 
porosity index corresponding to the effective porosity and 
that corresponding to the shale content. The latter can be 
estimated from a knowledge of the shale content (derived from 
the gamma ray curve) and a factor based on the porosity index 
for the case where the shale content approaches unity. Fig 2 
indicates this factor to be 55-60 per cent. 

In order to interpret the true resistivity in terms of fluid 
saturations, it is necessary to know how the shale is distributed 
in the sands, because shale laminations will have a different 
effect from disseminated shale particles. 

Accepting the relationship shown by the dashed line of Fig 
2, the resistivities of 100 per cent water-bearing shaly sands 
were calculated as a function of gamma ray deflections for 
sands with shale laminew (solid curve, Fig 3) and for sands 
with disseminated shales (dashed curve, Fig 3). Although 
the curves do not differ greatly, the dots representing the 
observations for the Third Burgan sand show better agreement 
with the solid curve, and hence the sand is inferred probably 
to have shale laminz. It is considered that in oil-boaring 
sands the relationship between gamma ray deflections and 


resistivity would not reveal so clearly the influence of the 
shales. 

In the absence of an induction log it would be very difficult 
to reach a conclusion regarding the shale distribution, because 
the R; values would not be very accurate. 
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TEMPERATURE NOTES ON THE BAHRAIN FIELD. W.H. Cotter. 


(The Bahrain Petroleum Co. Ltd.) 


When temperature data for wells completed in the Second 
Pay limestone of the Bahrain zones were corrected to a datum 
of 1900 ft subsea and plotted, the contours of equal tempera- 
ture broadly followed the structure contours, but gave a 
northerly shift of the temperature “crest” relative to the 


WELL No 63-5E(55) 


2010 


2020) 


2% pay 
Lmestone 2028 


TUBING T- iN | 
CASING FLOWING ON 4 14/64" BEAN 
TH P 425 PSI CHP 460 PY 
GAS RATIO 362 CUFT 


2050) 


2060) 


2070 


verrn FEET 


2090) 


LIMESTONE PERFORATIONS ~ DEPTH FEET. 


TUBING TAIL AT 2/27" 


2135» BoTTOW 2%! Pay LIMESTONE PERFORATIONS 


Fie 4 


JOURNAL OF THE INSTITUTE OF PETROLEUM 


| 
i 
Ri | | | | 
| R, | | 
| / 
2 3 4 5 b 
| | | 
| 
ae 2 | | | 
; 
| 
2040 
| 
: | | 
| \ 
| 
| 
eo 
220 
Py; 


SUMMARIES OF TECHNICAL NOTES PRESENTED AT THE KUWAIT SYMPOSIUM 73 


structural crest. From 300 ft subsea, where the temperature 
is fairly constant at 86° F, to the Bahrain zones, the tempera- 
ture gradient is about 3°3° F/100 ft. 

Wells in the north beyond two faults have water with a 
much higher sulphate content than those to the south in the 
same herizon. All waters in the lower pays of the Bahrain 
zones are rich in sulphate, suggesting that there may be move- 
ment of water via the faults to the northern sector. This may 
be of significance with regard to the temperature “ crest ”’ 
shift. 

From the Bahrain zones to the Arab zones the geothermal 
gradient is about 2° F/100 ft. 

Gas injection has been used for pressure maintenance in the 
Second Pay limestone of the Bahrain zones. A temperature 
increase of about 2° F has been observed near the topmost 
perforations (Fig 4) of those wells which have been reached by 
the gas, the arrival being shown by a sharp rise in gas/oil 
ratio. This temperature rise occurs after the gas reaches the 
well, and is attributed to free gas dissolving in the highly 
undersaturated crude. 
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A AnD Metruop or DETERMINING THE SALINITY 
or Some CrupvE Ors. P. Branson. (Qatar Petroleum 
Co. Ltd.) 


Oil with salt content in excess of 10 lb NaCl/100 bri tends 
to increase corrosion of refinery equipment. A simple test 
was, therefore, developed for determining the salt content of 
oil from Qatar. This involves shaking 75 ce of crude with 75 cc 
of distilled water containing three drops of Lissapol in a 
separating funnel for 10 minutes. After standing for 10 
minutes, 50 ce of water was drawn off and titrated with stand- 
ard silver nitrate solution. 

For salinities below 50 Ib NaCl/1000 bri the results were 
repeatable, and comparable with those obtained by Standard 
Method IP 77/51, although much lower in absolute accuracy. 

The method was not successful with Ain Zalah and Kirkuk 
crudes, because these gave emulsions from which only a few 
ce of clear water separated. A de-emulsifier other than 
Lissapol may be needed for the successful use of this method 
with these crudes. 
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GAS-LIQUID CHROMATOGRAPHY : EFFECT OF SUPPORT SIZE 
AND PROPORTION OF LIQUID PHASE ON COLUMN 
EFFICIENCY* 


By J. D. CHESHIRE and R. P. W. SCOTT + 


SUMMARY 


This paper describes the effect on efficiency of support size, proportion of liquid phase to support, and length 
ofcolumn, The conditions required for obtaining column efficiencies of 12,000 theoretical plates are given. The 
effect of size of charge on high-efficiency columns is also discussed. 


INTRODUCTION 


In their work on column efficiency Keulemans and 
Kwantes ! stated that the support size of 25 to 85 
B.S. mesh (particle diameter 599 to 178 u) produced 
columns of maximum efficiency. They also stated 
that the particle size affected eddy diffusion occurring 
in the column and expressed the eddy diffusion effect 
algebraically by the function, 2Ad,, originally derived 
by van Deemter, Zuiderweg, and Klinkenberg,? where 
2 is a characteristic of the packing and d, is the 
average particle diameter. This function contributes 
to the height of the theoretical plate and increases as 
the overall efficiency of the column decreases. The 
authors! explained that although it appears that a 
reduction of d, would increase the overall efficiency of 
the column, this is not always so. A point would be 
reached at which further reduction of d, would pro- 
duce irregular packing and channelling, resulting in 
an increase in the factor 4 and a corresponding fall in 
overall efficiency. 

The grade of 25 to 85 B.S. mesh, suggested by 
Keulemans and K wantes for column supports, has not 
been used by all workers. In their original work on 
gas-liquid chromatography, James and Martin * used 
Celite graded to an average particle diameter of 40 u. 
Adlard,‘ in his work on the polyethylene glycols, used 
supports graded to 85 to 120 B.S. mesh (178 to 124 y), 
and Mellor > stated broadly that the efficiency of a 
column increased as the particle diameter of the 
support was reduced. 

Widely differing quantities of liquid phase on the 
support have also been reported. Keulemans and 
Kwantes ! state that too high a proportion of liquid 
phase to support affects the mass transfer taking place 
in the column and reduces the overall efficiency. Too 
small a ratio of liquid phase, however, tends to intro- 
duce adsorption problems in connexion with the 
support itself. They have effected a compromise and 
state that the percentage of liquid phase on the 


support should be between 13 and 17 per cent w/w. 
Grant and Vaughan*® and Brooks, Murray, and 
Williams 7 used 23-1 per cent w/w liquid phase on 
support. James and Martin,’ Harrison,* and Adlard * 
all used 28-6 per cent w/w liquid phase on support, 
while Keppler, Dijkstra, and Schols ® in their work at 
high temperatures used between 9 and 33 per cent w/w 
liquid phase on the support. 

Thus an examination of the literature leads to the 
conclusions that there is still some uncertainty as to 
the best particle size to employ for the support and 
that there is little agreement on the quantity of liquid 
phase to be used on the support, values from 9 to 33 
per cent w/w having been reported. The work 
described in this paper was carried out to improve the 
separation of hydrocarbons. The correct particle size 
to be used for the support and the optimum proportion 
of liquid phase to support that is required to produce 
the maximum column efficiency were determined. 


PARTICLE SIZE OF SUPPORT 


Experimental 


Previous work had shown that with the particular 
method of packing used, Johns—Manville C22 insulat- 
ing brick was the most suitable material for the sup- 
port in gas-liquid chromatography. These bricks 
were ground, sieved to the following mesh sizes, and 
treated with 20 per cent Apiezon oil A as a liquid 
phase. 


Support grades Particle diameters 
30 to 60 B.S. mesh 500 to 251 pw 
60 to 100 251 to 
100 to 120 152 to 
120 to 160 124to 
100 to 200 152to 


The apparatus used was the routine equipment 
previously described by Scott,!° and the column was 
packed by vibration obtained from a low-speed elec- 


* MS received 15 November 1957. 


+ Benzole Producers Ltd. 
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tric motor. The column conditions used were as 
follows: 
Column length 5 ft 
Column diameter 3°6 mm 
Operating temperature 78°C 
Carrier gas . 25% nitrogen, 75% 
hydrogen 
Carrier-gas flow . 10 to 65 ml/min 
Column charge . 3 to 3-5 mg 


Chromatographs were obtained for n-heptane using a 
_ series of gas flows for each grade of support and the 
column efficiency calculated. The results are 
shown in Fig 1. 
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EFFICIENCIES FOR "-HEPTANE ON SUPPORTS OF VARIOUS 
PARTICLE SIZE 


Discussion 


Examination of Fig 1 ibis that the smaller the 
particle size of the support, the higher the resulting 
column efficiency at the optimum gas flow, which 
agrees with the results obtained by Mellor.° The 
optimum flow could not be reached for support graded 
120 to 160 B.S. mesh, as the pressure required to 
ensure the necessary flow of gas through the column 
was too high for the safe use of the glass injection 
system. It should be noted that the wide-range sup- 
port grade, 100 to 200 B.S. mesh, gives lower effi- 
ciencies than the narrow-range support grades, 100 to 
120 and 120 to 160 B.S. mesh. The results show that 
if the pressure at the inlet of the column is limited to 
a value that is safe for glass apparatus, the optimum 
size of support is the narrow grade of 100 to 120 B.S. 
mesh. 


QUANTITY OF LIQUID PHASE ON THE 
SUPPORT 


Experimental 
The apparatus was the same as that used in the 


investigation of optimum particle size. The column 
conditions were as follows: 
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Column length . 5ft 

Column diameter. 36mm 

Operating 78°C 

Carrier gas. 25% nitrogen, 75% 
hydrogen 

Carrier-gas flow 6 to 65 ml/min 

Support . Johns—-Manville C22 
insulating brick 

Grade of support 100 to 120 B.S. mesh 

Liquid phase . Apiezon oil A 

Column charge 0-2 to 3°5 mg 


A mixture of n-pentane, n-hexane, and n-heptane was 
chromatographed at various gas flows and the opti- 
mum flow, which gave the maximum column effi- 
ciency, was determined for columns carrying 2-5, 5, 
10, 15, 20, 25, and 30 per cent w/w of liquid phase on 
the support. A chromatogram of the following 
hydrocarbon mixture was then obtained for each 
proportion of liquid phase to support using the 
optimum flow previously determined: 


n-pentane 

n-hexane 

n-heptane 

n-octane 
methyleyclohexane 
benzene 

toluene 
ethylbenzene 
o-xylene 

p-xylene 
n-propylbenzene 
tsopropylbenzene 


A graph of the efficiencies obtained plotted against 
retention time for each proportion of liquid phase to 
support is shown in Fig 2. 

Discussion 

The results shown in Fig 2 indicate that, when 
employing a 5-ft column, 3-6 mm in diameter, to 
chromatograph substances that boil in the neighbour- 
hood of the column temperature, the optimum quan- 
tity of liquid phase to support that gives the maximum 
column efficiency is 20 per cent w/w. If the column 
is used at the optimum gas flow, this quantity of 
liquid phase will produce efficiencies of the order of 
2500 to 3500 theoretical plates. Such efficiencies will 
be obtained for substances with retention times lying 
between 5 and 60 minutes. 

The shapes of the curves in Fig 2 are also of con- 
siderable interest. Confining attention to curves for 
10 per cent or more of liquid phase to support, the 
efficiency increases rapidly with retention time during 
the first 5 minutes of the chromatogram. After this 
point the slope of the efficiency/retention time curve 
decreases suddenly, and, for a support carrying 30 
per cent w/w of liquid phase, the efficiency after 25 
minutes shows very small increase as the retention 
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time increases. The following equation, derived by 
van Deemter, Zuiderweg, and Klinkenberg,? and used 


EFFICIENCY IN THEORETICAL PLATES 
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Fie 2 
EFFECT OF PROPORTION OF LIQUID PHASE ON EFFICIENCY 


by Keulemans and Kwantes,! gives algebraic repre- 
sentation of all of the factors affecting column 
efficiency : 

H = 20d, + yD,/u + 8/n?.k’/(1 + k’)?.d?P/D,.u 

where H = height equivalent to a theoretical plate 
and, for a column of a given length, is 
inversely proportional to efficiency ; 
2 = a characteristic of the packing ; 

d, = average particle diameter ; 

y = a correction factor for the tortuosity of 
gas path; 

D, = diffusivity in gas phase; 

D, = diffusivity in liquid phase ; 

u = linear velocity of gas; 

k’ = K .F,/F,, where K = partition coefficient 
of solute and F; and F, = the fractions 
of cross section occupied by the liquid 
and gas phases respectively ; 

d,; = effective thickness of liquid film on the 
support. 


As mentioned before, the function 2d, accounts for 
the eddy diffusion occurring in the column. The 


function yD,/u accounts for the gaseous or molecular 
diffusion occurring in the column and 


8/n?. k/(1 +- k’)?. d?/D; 


accounts for the resistance to mass transfer. 

Under the conditions in which the results in Fig 2 
were obtained, Ad, and uw were constant and, for a 
given curve, D,, D;, and k'/(1 + k’)? were the only 
variables. The factor D, is associated with the 
gaseous diffusion effect only, and D; and k’/(1 + k’)? 
are associated with the resistance to mass transfer 
only. As the molecular weight, and thus the parti- 
tion coefficient and retention time, of the substances 
chromatographed increase, the value of D, decreases, 
resulting in an increase in the column efficiency. 
Conversely, as the molecular weight of the substance 
chromatographed increases, k’/(1 + k’)?.d;?/D, in- 
creases and results in a loss of efficiency. It appears 
from Fig 2 that at the beginning of the chromatogram 
the fact adversely affecting column efficiency is mainly 
gaseous diffusion, because increase in molecular weight 
of the substance chromatographed results in decrease 
in D, and the column efficiency increases with reten- 
tion time. However, a point is reached at which 
further increase in the molecular weight of the sub- 
stance results in the function k’/(1 + k’)? .d;?/D, 
becoming effective, which opposes the improvements 
resulting from the decrease in gaseous diffusion. It 
is concluded therefore that the change of slope in the 
efficiency/retention time curve occurs where the rate 
of change of resistance to mass transfer becomes of the 
same order as the rate of change of gaseous diffusion, 
and the rate of increase in column efficiency is thus 
reduced. The almost horizontal part of the curve for 
30 per cent w/w liquid phase on support indicates that 
theincrease in efficiency due to the reduction of gaseous 
diffusion is approximately equal to the increase in 
resistance due to mass transfer. However, this sud- 
den change of slope does not occur so early in the 
elution time curve, when 5 per cent w/w or less of 
liquid phase is present on the support. There is little 
or no change in the slope of the efficiency/retention 
time curve for elution times up to 12 minutes, and 
the efficiency continues to increase rapidly. This 
would be expected, as under these conditions d,? 
becomes so small that the resistance to mass transfer 
factor no longer becomes effective. Although sub- 
stances with elution times greater than 12 minutes 
could not be placed on the column due to low volatility 
at the operating temperature, it is obvious that this 
sudden change of slope would eventually occur if 
longer retention times could be realized. The effect 
of mass transfer has been considerably reduced, but, 
of course, cannot be completely eliminated. These 
results lead to two conclusions: 


(1) High-efficiency columns for chromato- 
graphing high-boiling liquids can best be pre- 
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pared by using short columns run at a low 
temperature with a low proportion of liquid 
phase to support. 

(2) High-efficiency columns for chromato- 
graphing low-boiling substances may be obtained 
by increasing the length of the column to give a 
longer chromatogram, but, at the same time, 
maintaining the same proportion of liquid phase 
to support, namely, 5 per cent w/w. 


CONDITIONS FOR HIGH EFFICIENCY 


Experimental 


To verify these conclusions, the following columns 
were constructed and their efficiencies determined: 


] 

Column 1 | Column2 | Column 3 
Length . 5 ft 10 ft 25 ft 
Diameter . 3-6 mm 3-6 mm 3-6 mm 
Operating temperature 78° 78°C 
Carrier gas 25% Ny/75% H, | 25% Ny/75% H, | 25% Ny/75% Hy 
Carrier-gas flow 37 ml min 37 ml/min 37 ml/min 
Support . | Johns—Manville | Johns—Manville | Johns—Manville 

| (22 insulating C22 insulating C22 insulating 

| brick brick brick 
Grade of support - | 100-120 B.S. 100-120 B.S. 60-100 B.S. 

| mesh mesh mesh 
Liquid phase. | Apiezon oil A Apiezon oil A Apiezon oil A 
Liquid phase to | 

support 5% wiw 5% w/w 5% wiw 

Column charge . 0-2-0-7 mg 0-2-0-8 mg 0-4 mg 
Inlet pressure . - | 20 Ib/sq. in. 22 Ib/sq. in. 35 Ib/sq. in. 


The support grade of 60 to 100 B.S. mesh was used 
for column 3 to reduce the back pressure at the 
column inlet. The injection system on column 3 was 
also modified to ensure safe use at the high inlet- 
pressure. 

These columns were used to chromatograph the 
following mixtures: 


T 
Mixture 1 | Mixture 2 Mixture 3 
| Boiling Boiling Boiling 
Substance | point | Substance point Substance point 

n-Pentane | 36-0 | n-Hexane 68-8 | n-Pentane . | 36-0 
n-Hexane | 68-8 | Benzene 80-1 | n-Hexane 68-8 
Benzene . 80-1 | n-Heptane . 98:4 | Methyleyelo- 
n-Heptane 98-4 | Methyleyclo- pentane 71-8 
Methyleyclo- | hexane 101-0 | Benzene 80-1 

hexane 101-0 | Toluene ° 110-8 | n-Heptane . 98-4 
Toluene . 110-8 Ethylbenzene | 136-2 | Methyleyelo- 
n-Octane. . | 125-6 | p-Xylene 138-4 hexane 101-0 
Ethylbenzene . 136-2 | o-Xylene | 144-0 | Toluene 110-8 
p-Xylene 138-4 | isoPropyl- | n-Octane . 125-6 
o-Xylene 144-0 benzene | 152-4 | Ethylbenzene | 136-2 
isoPropyl- n-Propyl- m-Xylene 137-0 

benzene 152-4 benzene . 158-6 | p-Xylene 138-4 
n-Propyl- Mesitylene . 164-6 | o-Xylene 144-0 

benzene 158-6 | p-Cymene 1773 
Mesitylene | 164-6 : 
n-Decane 174-0 | 
Hydrindene 178-0 | | | 
Indene . 182-0 | } } 
n-Undecane 195-8 | 


The optimum flow that gave the maximum effi- 
ciency was found to be approximately 37 ml/min for 
all three columns. Mixtures 1, 2, and 3 were chroma- 
tographed on columns I, 2, and 3 respectively. The 
chromatogram obtained from column 3 is shown in 
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Fig 3. The efficiency for each substance in all three 
chromatograms was calculated and the efficiencies, 
plotted against retention time, are shown in Fig 4. 
It should be stressed that the size of the column charge 
became very critical when using column 3. Previous 
work had shown that when using 20 per cent w/w 
liquid phase on the support, the column efficiency was 
not affected by the quantity of the charge, providing 
it was kept below 4 mg. However, when using 
columns | and 2 it was found that if the total charge 
exceeded 1 mg asymmetric peaks were formed and 
poor efficiency resulted. This effect was even more 
pronounced with column 3, for which a total charge in 
excess of 0-4 mg resulted in asymmetric peaks. This 
corresponded to a maximum of about 0-08 mg for any 
individual component present in the mixture. Due 
to the necessity for maintaining low column-charges, 
the sensitivity of the detector used in this work will 
be of interest and is stated as follows: 

A hydrocarbon produced 2 per cent of full-scale 
deflection of the recorder at a concentration of 0-05 ug 
per ml of carrier gas. This sensitivity was obtained 
with a long-term instability or drift of <1 per cent 
full-scale deflection and a short-term instability or 
noise of <0-3 per cent full-scale deflection. 


Discussion 


It may be seen from Fig 4 that a change of slope 
occurs in the efficiency/retention time curve obtained 
for all three columns, but that this inflection occurs at 


TABLE 1 
Efficiency in 
Substance theoretical 

plates 
n-Hexane 3620 
Benzene 3090 
n-Heptane . 4210 
Methyleyelohexane 4560 
Toluene 4340 
n-Octane 5440 
Ethylbenzene 4490 
p-Xylene 4370 
m- and p-Xylene . 3160 
o-Xylene 4980 
Cumene P 5250 
n- Propylbenzene ‘ 5360 
Mesitylene 5060 
p-Cymene 6520 


progressively higher efficiency values as the length of 
the column is increased. A 5-ft column run at 101° C 
may be used for chromatographing substances boiling 
up to 200° C with a column efficiency of up to 4000 
theoretical plates and a total elution time of 40 
minutes. A 25-ft column run at 78° C will produce 
chromatograms for substances boiling up to 144° C 
with an efficiency of about 12,000 theoretical plates 
and a total elution time of 120 minutes. Efficiencies 
of 8000 to 9000 plates were obtained from this column 
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for substances eluted in 10 minutes. No separation 
of m- and p-xylenes was shown on the 10-ft column. 
Some indication of the movement of the respective 
peaks of the meta and para isomers is given in Table I, 
which shows that the efficiency of the peak for the 
pure para isomer falls correctly in the series of effi- 
ciencies with respect to retention time. However, the 
peak for the mixture of the meta and para isomers 
shows a considerably lower efficiency than would be 
expected from their position in the retention series. 
This is due to the slight separation of the isomers. 
Separation of m- and p-xylene, however, is shown on a 
chromatogram obtained from the 25-ft column. It 
should be noted that the efficiency of the 25-ft column 
is approximately double that obtained from the 10-ft 
column and between three and four times that ob- 
tained from the 5-ft column. Wiebe ” states that 
column efficiency does not vary linearly with column 
length and that a point is reached at which increasing 
the length of the partition column results in little 
improvement of the overall efficiency. Wiebe, how- 
ever, was using 31 per cent w/w of liquid phase to 
support and was working under conditions in which 
mass transfer effect was reducing column efficiency 
considerably. It appears that if the quantity of 
liquid phase is reduced to 5 per cent w/w, the column 
efficiency varies almost linearly with the length of the 
column. A true linear relationship between effi- 
ciency and column length might possibly be achieved 
by using 3 or 4 per cent w/w liquid phase to support. 


CONCLUSIONS 


The particle size of the support that gives the 
highest column-efficiency commensurate with reason- 
ably low column inlet-pressures is the close-cut grade 
of 100 to 120 B.S. mesh. The adverse effect of liquid 
diffusion on column efficiency may be reduced by 


using 5 per cent w/w or less of liquid phase to support 
and, at the same time, increasing the column length 
to obtain practical retention values. A column of 
10,000 to 12,000 theoretical plates, suitable for the 
separation of hydrocarbons, may be obtained from a 
column 25 ft long, 3-6 mm in diameter, and packed 
with 5 per cent w/w of Apiezon oil on Johns—Manville 
C22 insulating brick graded to 60 to 100 B.S. mesh. 
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To obtain these high efficiencies, the total charge on 
the column must be less than 0-4 mg and a suitable 
detector must be employed with sufficiently high 
sensitivity and stability to produce the required 
chromatograms, 
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time increases. The following equation, derived by 
van Deemter, Zuiderweg, and Klinkenberg,? and used 


EFFICIENCY IN THEORETICAL PLATES 


° 5 10 15 20 25 30°) 
ELUTION TIME, min. 
Fie 2 
EFFECT OF PROPORTION OF LIQUID PHASE ON EFFICIENCY 


by Keulemans and Kwantes,! gives algebraic repre- 
sentation of all of the factors affecting column 
efficiency : 

H = 2d, + yD,/u + 8/n?.k'/(1 + k’)?.d?/D,.u 

where H = height equivalent to a theoretical plate 
and, for a column of a given length, is 
inversely proportional to efficiency ; 
2 = a characteristic of the packing; 

d, = average particle diameter ; 

y = a correction factor for the tortuosity of 
gas path; 

D, = diffusivity in gas phase; 

D, = diffusivity in liquid phase; 

u = linear velocity of gas; 

k’ = K . F,/F,, where K = partition coefficient 
of solute and F; and F, = the fractions 
of cross section occupied by the liquid 
and gas phases respectively ; 

d; = effective thickness of liquid film on the 
support. 


As mentioned before, the function 2ad, accounts for 
the eddy diffusion occurring in the column. The 


function yD,/u accounts for the gaseous or molecular 
diffusion occurring in the column and 


8/n?. k/(1 +k’). d,?/D; 


accounts for the resistance to mass transfer. 

Under the conditions in which the results in Fig 2 
were obtained, Ad, and u were constant and, for a 
given curve, D,, D;, and k’/(1 + k’)? were the only 
variables. The factor D, is associated with the 
gaseous diffusion effect only, and D; and k’/(1 + k’)? 
are associated with the resistance to mass transfer 
only. As the molecular weight, and thus the parti- 
tion coefficient and retention time, of the substances 
chromatographed increase, the value of D, decreases, 
resulting in an increase in the column efficiency. 
Conversely, as the molecular weight of the substance 
chromatographed increases, k’/(1 + k’)®.d,;?/D, in- 
creases and results in a loss of efficiency. It appears 
from Fig 2 that at the beginning of the chromatogram 
the fact adversely affecting column efficiency is mainly 
gaseous diffusion, because increase in molecular weight 
of the substance chromatographed results in decrease 
in D, and the column efficiency increases with reten- 
tion time. However, a point is reached at which 
further increase in the molecular weight of the sub- 
stance results in the function k’/(1 + k’)? .d;?/D, 
becoming effective, which opposes the improvements 
resulting from the decrease in gaseous diffusion. It 
is concluded therefore that the change of slope in the 
efficiency/retention time curve occurs where the rate 
of change of resistance to mass transfer becomes of the 
same order as the rate of change of gaseous diffusion, 
and the rate of increase in column efficiency is thus 
reduced. The almost horizontal part of the curve for 
30 per cent w/w liquid phase on support indicates that 
theincrease in efficiency due to the reduction of gaseous 
diffusion is approximately equal to the increase in 
resistance due to mass transfer. However, this sud- 
den change of slope does not occur so early in the 
elution time curve, when 5 per cent w/w or less of 
liquid phase is present on the support. There is little 
or no change in the slope of the efficiency/retention 
time curve for elution times up to 12 minutes, and 
the efficiency continues to increase rapidly. This 
would be expected, as under these conditions d,? 
becomes so small that the resistance to mass transfer 
factor no longer becomes effective. Although sub- 
stances with elution times greater than 12 minutes 
could not be placed on the column due to low volatility 
at the operating temperature, it is obvious that this 
sudden change of slope would eventually occur if 
longer retention times could be realized. The effect 
of mass transfer has been considerably reduced, but, 
of course, cannot be completely eliminated. These 
results lead to two conclusions: 


(1) High-efficiency columns for chromato- 
graphing high-boiling liquids can best be pre- 
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pared by using short columns run at a low 
temperature with a low proportion of liquid 
phase to support. 

(2) High-efficiency columns for chromato- 
graphing low-boiling substances may be obtained 
by increasing the length of the column to give a 
longer chromatogram, but, at the same time, 
maintaining the same proportion of liquid phase 
to support, namely, 5 per cent w/w. 


CONDITIONS FOR HIGH EFFICIENCY 
Experimental 


To verify these conclusions, the following columns 
were constructed and their efficiencies determined: 


| Column 1 | Column 2 Column 3 
Length - 5 ft 10 ft 25 ft 
Diameter . 3-6 mm 3-6 mm 3-6 mm 
Operating temperature | 101°C 78° 0 78°C 
Carrier gas 4 . | 25% N4/75% Hy | 25% Ny/75% H, | 25% N,/75% Hy 
Carrier-gas flow 37 ml/min 37 ml/min 37 ml/min 
Support . . | Johns—Manville | Johns~Manville | Johns—Manville 

| O©22 insulating O22 insulating C22 insulating 

brick brick brick 

Grade of support |, 100-120 B.S. 100-120 B.S. 60-100 B.S. 

| mesh mesh mesh 
Liquid phase. . | Apiezon oil A Apiezon oi) A Apiezon oil A 
% Liquid phase to | 

support 5% wiw 5% w/w 5% 

Column charge . - | 0-2-0-7 mg 0-2-0-8 mg 0-4 mg 
Inlet pressure . 20 Ib/sq. in. 22 lb/sq. in. 35 Ib/sq. in. 


The support grade of 60 to 100 B.S. mesh was used 
for column 3 to reduce the back pressure at the 
column inlet. The injection system on column 3 was 
also modified to ensure safe use at the high inlet- 
pressure. 

These columns were used to chromatograph the 
following mixtures: 


Mixture 1 | Mixture 2 Mixture 3 
| Boiling | Boiling Boiling 
Substance | point | Substance point Substance point 
| °C 
n-Pentane | 36-0 | n-Hexane . 68-8 | m-Pentane . 36-0 
n-Hexane | 68-8 | Benzene : 80-1 | n-Hexane. 68-8 
Benzene . | 80-1 | n-Heptane . 98-4 | Methyleyclo- 
n-Heptane | 98-4 | Methyleyelo- pentane . 71:8 
Methylcyclo- | hexane . 101-0 | Benzene 80-1 
hexane ’ 101-0 Toluene > 110-8 | n-Heptane . 98-4 
Toluene . 110-8 Ethylbenzene | 136-2 | Methyleyclo- 
n-Octane. 9 125-6 | p-Xylene . | 138-4 hexane .| 101-0 
Ethylbenzene . 136-2 | o-Xylene - | 144-0 | Toluene . 110-8 
p-Xylene 138-4 | isoPropyl- n-Octane . 125-6 
o-Xylene . 144-0 benzene . 152-4 | Ethylbenzene | 136-2 
isoPropyl- n-Propyl- m-Xylene . 137-0 
benzene | 152-4 benzene. 158-6 | p-Xylene . 138-4 
n-Propyl- | Mesitylene . 164-6 | o-Xylene. 144-0 
benzene | 158-6 | p-Cymene ./| 
Mesitylene | 164-6 
n-Decane 174-0 
Hydrindene | 178-0 
Indene . 182-0 
n-Undecane 195-8 


The optimum flow that gave the maximum effi- 
ciency was found to be approximately 37 ml/min for 
all three columns. Mixtures 1, 2, and 3 were chroma- 
tographed on columns I, 2, and 3 respectively. The 
chromatogram obtained from column 3 is shown in 
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Fig 3. The efficiency for each substance in all three 
chromatograms was calculated and the efficiencies, 
plotted against retention time, are shown in Fig 4. 
It should be stressed that the size of the column charge 
became very critical when using column 3. Previous 
work had shown that when using 20 per cent w/w 
liquid phase on the support, the column efficiency was 
not affected by the quantity of the charge, providing 
it was kept below 4 mg. However, when using 
columns | and 2 it was found that if the total charge 
exceeded 1 mg asymmetric peaks were formed and 
poor efficiency resulted. This effect was even more 
pronounced with column 3, for which a total charge in 
excess of 0-4 mg resulted in asymmetric peaks. This 
corresponded to a maximum of about 0-08 mg for any 
individual component present in the mixture. Due 
to the necessity for maintaining low column-charges, 
the sensitivity of the detector used in this work will 
be of interest and is stated as follows: 

A hydrocarbon produced 2 per cent of full-scale 
deflection of the recorder at a concentration of 0-05 ug 
per ml of carrier gas. This sensitivity was obtained 
with a long-term instability or drift of <1 per cent 
full-scale deflection and @ short-term instability or 
noise of <0-3 per cent full-scale deflection. 


Discussion 
It may be seen from Fig 4 that a change of slope 


occurs in the efficiency/retention time curve obtained 
for all three columns, but that this inflection occurs at 


TABLE 
Efficiency in 
Substance theoretical 

plates 
n-Hexane. 3620 
Benzene ‘ 3090 
n-Heptane . 4210 
Methylcyclohexane ‘ 4560 
m- and p-Xylene . 3160 
Cumene ‘ 5250 
n-Propylbenzene . 4 5360 
Mesitylene . 5060 
p-Cymene . 6520 


progressively higher efficiency values as the length of 
the column is increased. A 5-ft column run at 101° C 
may be used for chromatographing substances boiling 
up to 200° C with a column efficiency of up to 4000 
theoretical plates and a total elution time of 40 
minutes. A 25-ft column run at 78° C will produce 
chromatograms for substances boiling up to 144° C 
with an efficiency of about 12,000 theoretical plates 
and a total elution time of 120 minutes. Efficiencies 
of 8000 to 9000 plates were obtained from this column 
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for substances eluted in 10 minutes. No separation 
of m- and p-xylenes was shown on the 10-ft column. 
Some indication of the movement of the respective 
peaks of the meta and para isomers is given in Table I, 
which shows that the efficiency of the peak for the 
pure para isomer falls correctly in the series of effi- 
ciencies with respect to retention time. However, the 
peak for the mixture of the meta and para isomers 
shows a considerably lower efficiency than would be 
expected from their position in the retention series. 
This is due to the slight separation of the isomers. 
Separation of m- and p-xylene, however, is shown on a 
chromatogram obtained from the 25-ft column. It 
should be noted that the efficiency of the 25-ft column 
is approximately double that obtained from the 10-ft 
column and between three and four times that ob- 
tained from the 5-ft column. Wiebe ™ states that 
column efficiency does not vary linearly with column 
length and that a point is reached at which increasing 
the length of the partition column results in little 
improvement of the overall efficiency. Wiebe, how- 
ever, was using 31 per cent w/w of liquid phase to 
support and was working under conditions in which 
mass transfer effect was reducing column efficiency 
considerably. It appears that if the quantity of 
liquid phase is reduced to 5 per cent w/w, the column 
efficiency varies almost linearly with the length of the 
column. A true linear relationship between effi- 
ciency and column length might possibly be achieved 
by using 3 or 4 per cent w/w liquid phase to support. 


CONCLUSIONS 


The particle size of the support that gives the 
highest column-efficiency commensurate with reason- 
ably low column inlet-pressures is the close-cut grade 
of 100 to 120 B.S. mesh. The adverse effect of liquid 
diffusion on column efficiency may be reduced by 


using 5 per cent w/w or less of liquid phase to support 
and, at the same time, increasing the column length 
to obtain practical retention values. A column of 
10,000 to 12,000 theoretical plates, suitable for the 
separation of hydrocarbons, may be obtained from a 
column 25 ft long, 3-6 mm in diameter, and packed 
with 5 per cent w/w of Apiezon oil on Johns—Manville 
C22 insulating brick graded to 60 to 100 B.S. mesh. 
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To obtain these high efficiencies, the total charge on 
the column must be less than 0-4 mg and a suitable 
detector must be employed with sufficiently high 
sensitivity and stability to produce the required 
chromatograms. 
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Barclays Bank Building « 73 Cheapside «+ London, England 
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TOOL COMPANY. 
WORLD STANDARD 


this Filter talk will 
drive round the bend! 


Sorry, old fruit. I was only trying to be helpful. 


Helpful! Here am | practising for the only real sporting event of the year and 
all you can talk about is some wretched little gadget on one of your 
blasted hairyplanes. Can't you see | need to concentrate? 


It’s not a little gadget and they'd never put it 
on a plane—it’s bigger than this Renault. 
They put it in the line between the storage 
and aircraft tanks to make sure 

the plane gets clean, water-free fuel. 


That's the trouble with your infantile 
modern engines-—they have to be fed on 
strained foods or they get indigestion. 
Now look at Dotty... 


Exactly. Do you really think she 

enjoys having that steam pouring 
through her exhaust? Now if you 

took a Simmonds Fram Separator Filter 
in tow you could have 100%, water 
separation and complete filtration 

of solids down to 5 microns. 


It's tempting, very 
tempting ... but it 
wouldn't be fair on 
the other competitors. 
Now d'you mind if we 
stop a minute... 

| think | really will 
have to get the 

rest of that 

bird’s nest out 

of the tank. 


the point of the argument... FRAM SEPARATOR FILTERS 
for clean water-free fuel 


For more details of this and other Simmonds products 
SEND FOR LITERATURE to 


SIMMONDS AEROCESSORIES LTD 
Byron House, 7-8-9 St. James’s Street, London, $.W.1. 


Head Office and Works: Treforest, Pontypridd, Glamorgan 


Also Birmingham, Glasgow, Stockholm, Copenhagen, Ballarat, 
Sydney, Johannesburg, Amsterdam and New York 


A MEMBER OF THE FIRTH CLEVELAND GROUP (FG) 
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Two famous companies join together 


YORKSHIRE IMPERIAL METALS LIMITED 


YORKSHIRE 


IMPERIAL 


.... the Best of both 


Yorkshire Imperial Metals Ltd. is the 
name of the new Company which combines 
the copper and alloy tube, fittings and 
plate activities of The Yorkshire Copper 
Works Ltd. and Imperial Chemical In- 
dustries Ltd., and also the plastic tube and 
plastic fittings interests of The Yorkshire 
Copper Works Ltd. 


The new Company thus draws upon 
unequalled research, development, technical 
and production services. The deployment 
of these in a single enterprise will promote 
economy, efficiency and_ service to 
customers, both large and small, in this 
country and abroad. 


The new Company means to keeps its 
friends, and to make new ones, by offering 
an unrivalled service in the supply of high 
quality tubes, plates and fittings, designed 
and manufactured to suit individual service 
conditions. 


The emphasis on service is all-important 
—the new Company intends to remain 
personal in its customer relationship. 


This relationship, built up over the years 
on terms of friendly co-operation, personal 
contact and real service, is to be cherished. 
No problem is ever too small to merit the 
Company’s undivided attention—it is here 
to help, and it is willing and able to help 
more than ever before. 


Whatever your requirements in copper 
and alloy tubes, plates or fittings, Yorkshire 
Imperial wish to be of service to you in 
many ways. Whatever the job you want 
them to do, Yorkshire Imperial offer pro- 
ducts of a quality which cannot be surpassed 
anywhere in the world, personal attention 
to your requirements, fair prices and delivery 
on time. 


specialising in the design and manufacture 
of non-ferrous tubes, plates and fittings— 
and plastic tubes and fittings 


HEAD OFFICE—LEEDS 


Factories ;— 


YORKSHIRE 


COPPE 
LEEDS 


. KIRKBY WORKS, . ALLEN EVERITT 


R WORKS, LIVERPOOL 


Staffs. 


. FYFFE WORKS, 
WORKS, SMETHWICK, CAROLINA PORT, 
DUNDEE 


ANSON UNITS, . 
CASTLEFORD, 
Yorks. 


. LANDORE WORKS, . 


SWANSEA GLASGOW 
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PRESERVING A WAY OF LIFE 


AN OIL PIPELINE cannot be considered solely as a notable engineering accomplishment 
although it is a tribute to the wide range of skills of a modern engineer, presenting 
complex economic, technical and geographical problems. More than that—it is a 
symbol of our Western civilization, and a promise to generations to come of its 
continuation and development. 
A pipeline preserves our Way of Life—feeding Industry, speeding transport, lubri- 
cating, motivating and making mass production possible. It means insecticides, 
fertilizers, fibres, food and life-saving drugs. The pipeline is vital in terms of trade, 
vital in terms of humanity. 
Pipelines themselves must be preserved—preserved against the ravages of corrosion. 
and for this purpose engineers have chosen ** BITUMASTIC ” coatings the world 
, over and notably in the Middle East. We are proud of their reliability, their high 
uniform standard of quality, and in our minor role, of the part we play in preserving 
the life-lines of those who will follow us. 
Many leading Oil Companies specify *“*‘ BITUMASTIC” coatings for protecting 
pipelines. Send for a copy of our Pipeline Data Booklet for full particulars of 
** BITUMASTIC ” products and their use in the Oil Industry. 


SPECIALISTS IN ANTI-CORROSIVE COATINGS 
SPECIALISTS IN PIPELINE PROTECTION 


H 134 WAILES DOVE BITUMASTIC LTD + HEBBURN + CO. DURHAM 
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Above 
Suspicion 


These best quality Mild Steel 


corrosion-resisting Tanks 


under 


are made in capacities 


from 250 gallons to 5,000 gallons. 


ground 


Larger sizes to order. 


A series of small capacity 
tanks ready for despatch 


UNDERGROUND in 


STORAGE TANKS Mild 


Available from stock Steel 


made by W. P. Butterfield Ltd P.O. Box 38 Shipley, Yorks Tel. 52244 (8 lines) 

Branches LONDON Tel. HOLborn 2455 (4 lines) BIRMINGHAM Tel. EAS 0871 & EAS 2241 BRISTOL Tel. 27905 LIVERPOOL 
Tel. Central 0829 MANCHESTER Tel. Blackfriars 9417 NEWCASTLE-ON-TYNE Tel. 23823 GLASGOW Tel. Central 7696 
BELFAST N.!. Tel. 57343 DUBLIN Tel. 77232 
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TO MAINTAIN CLOSE PERSONAL CONTACT 
WITH OUR CUSTOMERS, DRESSER OPERATES WORLD WIDE 


DRESSER MANUFACTURING COMPANY, LIMITED 
1211 Bathurst Street 
Toronto 4, Ontario, Canada 


CLARK BROS. CO. 
Ling House, Dominion Street 
London, E.C. 2, England 


PETRO-TECH SERVICE, C.A. 
Apartado 2728 
Caracas, Venezuela 


CLARK COMPRESSOR CO. LTD. 
1210 11th Avenue, W. 
Calgary, Alberta, Canada 


DRESSER MANUFACTURES (ENGLAND) LIMITED 
39 Victoria Street 
London, S.W. 1, England 


ROOTS-CONNERSVILLE BLOWER 
(CANADA), LTD. 
629 Adelaide St. West 
i Toronto 3, Ontario, Canada ; 
CLARK-ITALY S.p.A. | LANE-WELLS CANADIAN COMPANY 
Via Tembien 41, Rome, Italy | Room 14, 10548 82nd Avenue SECURITY ENGINEERING CANADA LTD. 
P.O. Box 4267 
! 
! 


Edmonton, Alberta, Canada 
DRESSER A.G. Edmonton, Alberta, Canada 


ce Mihlebachstrasse 43, Ziirich, Switzerland MAGNET COVE BARIUM CORPORATION, LTD. 
525-A Seventh Avenue, West 
Calgary, Alberta, Canada 


SECURITY ROCK BITS LTD.* E 
Whitworth Street 
Openshaw, Manchester 11, England 


DRESSER LIMITED, C.A. 
Apartado 2728 


MAGCOBAR DE MEXICO, S.A. 
Apartado 1742, Monterrey, N.L., México SOCIETE FRANCAISE 
DES INDUSTRIES DRESSER, S.A. 


11, rue Auber, Paris, France 


DRESSER INTERNATIONAL, INC.: 
Viamonte 867, Buenos Aires, Argentina 
Edif. Radio Continente, Ave. México 
Los Caobos, Caracas, Venezuela 


MAGCOBAR DE VENEZUELA, C.A. 
Edif. Luz Eléctrica, Ave. Urdaneta 
Caracas, Venezuela SOUTHWESTERN INDUSTRIAL 


Clark Bros. Pan American Div.— 
ELECTRONICS (CANADA) LTD. 


hee, Paseo de la Reforma 95-1102; 


Southwestern Industrial Electronics Div.—Edison 40-1; _ PACIFIC PUMPS OF CANADA, LIMITED 434 Seventh Ave., East 
Ideco Div.—813 Edificio Internacional, | 9707 63d Avenue Calgary, Alberta, Canada 
Paseo de la Reforma 1, México D. F., México | Edmonton, Alberta, Canada 


*Majority owned 


Keyed to your needs 
and linked to give 
DRESSER Plus #service 


“One of the Dresser Industries” is a well known phrase identifying 
the family relationship of 12 great Dresser companies which 
supply equipment and technical services to the oil, gas, 

chemical, electronic and other industries. 


Each Dresser company is independently managed for maximum 
efficiency. Each possesses global experience, engineering talent and 
modern manufacturing facilities keyed to the highly technical 
needs of the industries it serves. Linked by common ownership, 
Dresser companies can be quickly mobilized to meet special 
conditions calling for specialized research, engineering or 
production facilities. This coordinated teamwork is known as 
Dresser plus gservice. It’s an extra value provided promptly by 
men with imagination. The Dresser plus gis your 

positive assurance of superior equipment and technical services — 
the standard of comparison the world over. 


Tomorrow’s progress planned today 


ASTRIES, INC. 


+ GAS 
EQUIPMENT AND | cHEmicaL 


TECHNICAL SERVICES | ELECTRONIC 
INDUSTRIAL 


REPUBLIC NATIONAL BANK BUILDING, DALLAS 21, TEXAS 


These are the Dresser Industries: 


CLARK BROS. CO.—compressors DRESSER-IDECO COMPANY-—stcel structures DRESSER MANUFACTURING DIVISION 
—couplings THE GUIBERSON CORPORATION — oil tools IDECO, INC.—drilling rigs LANE-WELLS COMPANY -— teciinical 
oilfield services © MAGNET COVE BARIUM CORPORATION — drilling mud ¢ PACIFIC PUMPS, INC.—pumps « ROOTS-CONNERS- 
VILLE BLOWER DIVISION—blowers ©¢ SECURITY ENGINEERING DIVISION—drilling bits © SOUTHWESTERN INDUSTRIAL 
ELECTRONICS —electronic instrumentation WELL SURVEYS, INC.—nuclear and electronic research and development. 
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Rocksil 


resilient 


too 


ROCKSIL — for either hot or cold insulation — 
has many important advantages to offer. 
Probably the most outstanding is its 
combination of strength with resilience —‘rigid’ 
resilience that absorbs impact and stands up well 
to considerable deflection of the surface it is 
protecting. In bonded slab form of various 
thicknesses and densities this makes Rocksil 
ideal for casing and structural insulation, 

It will not settle under continuous vibration 
even as loose infill. Rocksil is produced from 
natural Scottish rock and is odourless, 
chemically inert and sulphur-free, 


non-hygroscopic and rot-proof. 


ocksil 


Full information available from the manufacturers : 
THE CAPE ASBESTOS COMPANY LTD. vv 
114 & 116 Park Street, London W. 1. Tel: Grosvenor 6022 “ey 


ROCK WOOL INSULATION 


Glasgow: Eagle Buildings, 217 Bothwell Street, Glasgow C.2. 
Telephone : Central 2175. 


Manchester: Floor D, National Buildings, St. Mary's Parsonage, 
Manchester 3. Telephone : Deansgate 6016-7-8. 
Birmingham: || Waterloo Street, Birmingham 2. 
Telephone : Midland 6565-6-7, 


Newcastle: 19 & 20 Exchange Buildings, Newcastle-upon-Tyne. 
Telephone : Newcastle 20488, 
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BEYOND 
ALL QUESTION 


TIGHT 


NEWMAN -McEVOY VALVES GIVE THE HIGHEST EFFICIENCY 


Newman-McEvoy Valves are absolutely tight 
at all times, under any and every condition. 
Automatic continuous self-sealing, self- 
aligning gates and  non-distorting seats 
further provide perfect fit and tightness. 
Newman-McEvoy Valves are very easy to 


operate and are thoroughly economical in 

maintenance. It is these prime qualities of the 

valves, backed by other exceptional features, 

that account for their presence on the wellhead 

equipment of some of the greatest oilfields. 
Fullest details and literature on request. 


Newman, Hender & Co. Ltd 


CODCHESTER STROUD, GLOS. 
Telephone : Lw Telegrams : VALVES, STROUD 
n 


Telex 43-220 
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For present AND FUTURE A\S.T.M. and LP. 


Viscosity Specifications 


THE NEW X.27 MAJOR 
VERY CONSTANT 
TEMPERATURE BATH 
ACCURACY +0-010°C AT 100°C 


(and even better at lower temperatures) 


Takes 8 viscometers in prelevelled lid. 
The automatic boost heating allows quick rise - 
to any preset temperature and then automatic 
stabilisation. 


TOWNSON AND MERCER LIMITED 


CROYDON ENGLAND 
Telephone: THOrnton Heath 6262 


We Specialise in 


ALL TYPES OF 


STEEL 
STRUCTURES 


for the 


OIL INDUSTRY 


LAN 


Ay, 


‘KELVIN’ iron and ‘MAINSTEEL’ PALISADING, 
RAILINGS, GATES 


LONDON OFFICE WORKS AND REGISTERED OFFICE 
VINCENT HOUSE, VINCENT SQUARE, S.W.1 CLYDESDALE IRON WORKS, POSSILPARK 
GLASGOW, N.2 


Telephone: Victoria 8375 Telegrams: Kelvin Sowest, London Telephone: Possil, 8381 Telegrams: Kelvin, Glasgow 


CALCUTTA CHITTAGONG _ NATROBI 


xiv 


: 
= 


This UOP Platformer 
serves the German petro- 
leum market with Plat- 
formate (90 F-1 octane) 
at a rate of 3,000 barrels 
per stream day. 


UNIVERSAL OIL 


30 Algonquin Road, 

Des Piaines, lilinois, U.S.A. 
Representative in England: F. A. Trim, 
Bush House, Aldwych, London, W. C. 2 
More Than Forty Years Of Leadership 
in Petroleum Refining Technology 


*Registered Trademarks Of The 
Universal Oil Products Company 


* PRODUCTS COMPANY 


PLATFORMING* 
PROCESS 


xv 


ae Like so many other refiners the world over oe : 
ss at Misburg, Germany serves the needs of its customers with the oa: 


FLOW BOX 


PRESSURE 0-200'w.g. 
8 TEMPERATURE 50°-150°r 


SET POINT 
VALVE POSITION 


4-INCH CHART (same recording width as 
a 12-inch circular chart). Instrument cut-out 
size only 63” » 63”. 


@ 4 PRESSURE MOVEMENTS (2 Recording 
plus 2 Indicating). Ideal for Automatic 
Control. 


@ automatic CHART PICK-UP and 
ECONOMISER, (Patented) Chart lasts for 
1,000 hours at }” per hour. 


@ TRUE RECTILINEAR MOVEMENT 
of pens and Indicators. No more “arc” 
traces. 


An informative Leaflet No. S91 cover- 
ing the comprehensive range o 
Electrofio Recorders and Recording 
Control Stations is available on 
request. 


ELECTROFLO METERS COMPANY LIMITED Head Office: Abbey Road, Park Royal, London, N.W.10 
Factories; Abbey Road and Chase Estate, Park Royal and Maryport, Cumberland Telephone: Elgar 7641/8. 
Grams & Cables; Elflometa, London, Telex. Teler: No, 2-3196 


RECORDER Made by N ATOR 
: 
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FROM NATURAL GAS OR NAPHTHA 


Through special arrangement 
with Societe Belge de Azote 
(SBA), The Kellogg Inter 
national Corporation now can 
provide a continuous acetylene 
production process using naph- 
tha or natural gas as feedstock. 
This improved process pro- 
mises economies and_ results 
previously not obtainable. 


Revolutionary Burner 
The heart of this process is the 
SBA-developed burner. Easy to 
control, and requiring only 
minimum attention, it is stable 
and operates continuously — 
producing a constant yield. The 
formation of coke and_ tarry 
materials in the burner is virtu- 
ally non-existent. Any trace of 
these materials is automatically 
removed by special separation 
equipment. 


Kellogg Cooperation 
Through cooperative research 


and development by SBA and 
Kellogg, considerable engineer- 
ing information has been de- 
veloped on both the pyrolysis 
and purification systems. The 
results of this joint effort, 
available only from Kellogg, 
assure accurate and thoroughly 
reliable plant designs. 


Process Highlights 
Using naphtha feedstock, a 
wide range of ethylene to acety- 
lene ratios may be achieved. 
This permits simultaneous pro- 
duction of ethylene and acety- 
lene with great flexibility. In 
the case of natural gas feed, 
acetylene is the primary pro- 
duct—with negligible quanti- 
ties of heavy impurities. If de- 
sired, the plant can be designed 
to operate on either feedstock. 


Acetylene Purification 
In the process train, the heavy 


acetylenes and hydrocarbons 
are removed first, leaving a 
stream containing. acetylene 
and lighter materials. By using 
an SBA-developed ammonia 
absorption system, acetylene 
and ethylene are recovered from 
hydrogen and the other products 
of combustion, which in them- 
selvesare valuable by-products. 
The acetylene, 99.5°,, pure or 
higher, andethylenearestripped 
from the ammonia. Ammonia 
makes this process safer than 
others because of its depressing 
effect on the explosibility of 
acetylene. 


For More iIntormation 


Chemical companies and petro- 
leum refiners interested in utilis- 
ing this advanced acetylene 
production process are invited 
to discuss their specific needs 
in person with The Kellogg 
International Corporation. 


Kellogg International Corporation 


KELLOGG HOUSE - 17-10 CHANDOS STREET, CAVENDISH SQ., 


SOCIETE KELLOGG PARIS 


THE CANADIAN KELLOGG COMPANY LTD - TORONTO 
KELLOGG PAN AMERICAN CORPORATION = NEW YORK 
COMPANHIA KELLOGG BRASILEIRA * RIO DE JANEIRO 
COMPANIA KELLOGG DE VENEZUELA © CARACAS 


Subsidiaries of THE M. W. KELLOGG COMPANY 
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BROTHERHOOD 


VERTICAL AND HORIZONTAL 


STEAM TURBINES 


FOR DRIVING PUMPS, ETC. 


Wide range—All types. 
ee | Over 50 years’ experience. 

Hundreds in hand— 

thousands in service. 


BROTHERHOOD 


COMPRESSORS 


Air, GasandRefrigerating. 


The widest range in the 
British Empire—made tosuit 
your requirements. 
Thousands in service. 


BROTHERHOOD 
GENERATING SETS 


Turbine driven up to 

11,000 kw. 

Engine driven up to 340 kw. 
Scores in hand, hundreds 
in service. 


BROTHERHOOD 
REFRIGERATING PLANT 


Ammonia, CO,, Freon, SO,, 
Methy! Chloride. Wide range 
—single and double acting— 


one or more stages. 


Made to measure for 
special duties. 


also Manufacturers of all kinds of 
PLANT TO CUSTOMERS’ OWN DESIGNS | 


WHY NOT SEND YOUR PROBLEMS TO US? 


We shall be pleased to i gate them confidentially 
without commitment 


“COMPRESSOR & POWER PLANT SPECIALISTS FOR NEARLY A CENTURY 


STANDARD METHODS 


FOR 


TESTING PETROLEUM 


AND 


ITS PRODUCTS 


(Excluding Engine Test Methods for Rating Fuels) 


(SEVENTEENTH EDITION, 1958) 


788 pages Illustrated 


Price 40s. post free 


Obtainable from 


The Institute of Petroleum 
61 New Cavendish Street, 
London, W.1 
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and associated equipment 


ELECTRIC MOTORS 


FOR THE OIL INDUSTRY 


Metrovick Engineers have made an extensive 
study of motor, control gear and switchgear 
problems in the Oil Industry. Results of this 
experience are crystallised and arranged for easy 
reference in this fully illustrated 197 page book. 
Information given includes a selected choice of 
motors and associated equipment both flame- 
proof and non-flame proof, with rating lists and 
dimension tables. Useful general details of motor 
drives and maintenance are also included. 
Senior personnel associated with oil and chemi- 
cal industries are invited to write for publication 
7765/2. 


METROPOLITAN -VICKERS 


ELECTRICAL CO LTD. TRAFFORD PARK MANCHESTER 
‘ 


An A.E.1. Company 


Motors for the Oil Industry » 
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45 cubic metre Pressure Vessel 


made for I.C.I. Billingham 
Division. Length 30 feet. 
Diameter 10 feet. Fabricated 
from Kynal M.35 1. Alu- 
minium Alloy 3” thick. One 
end dished and a rolled body 
cone and skirt. Argon Arc : 
welded throughout. — 


MARSTON EXCELSIOR in the Service of Industry 


Marston’s have unrivalled experience in the fabrication of light alloys—of all shapes and sizes. 
Their products have earned a reputation for efficiency and reliability that is world-wide. 


* Light Alloy Fabrication 


* Specialised Engineering Assemblies 


* Laminated Plastic Components 
* Flexible Tanks 
* Radiators and Heat Exchangers 


Tubular Heat Exchangers ba Oil Tank 


5 pass cross contra type made for 
I.C.1. Billingham Division. Length 
184 feet. Diam. 235” inside. Header 
tanks are Aluminium Castings 
of BS.1490 LM.20 alloy 3” thick. 
Body rolled and welded from 
PA.19 Aluminium Alloy Plate }” 
thick. 

Tested to 7§ p.s.i. 


of 2500 gallon capacity, made 
for I.C.I. Paints Division. 
Length 14 feet. Diameter 7 feet. 
Fabricated from Kynal_ P3. 
Aluminium alloy }” thick with 
dished ends and a rolled plate 
body. Argon Arc welded 
throughout. Tested to withstand 
an 18 foot head of water. 


MARSTON EXCELSIOR LIMITED 


(A subsidiary company of Imperial Chemical Industries Ltd.) 
FORDHOUSES, WOLVERHAMPTON. Tel: Fordhouses 3361 
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**Newallastic’’ bolts and studs have qualities which are absolutely unique. 


They have been tested by every known device, and have been proved to 


be stronger and more resistant to fatigue than bolts or studs made by 


the usual method. 
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45 cubic metre Pressure Vessel 


made for I.C.I. Billingham 
Division. Length 30 feet. 
Diameter 10 feet. Fabricated 
from Kynal M.35 1. Alu- 
minium Alloy }” thick. One 
end dished and a rolled body 
cone and skirt. Argon Arc 
welded throughout. 


MARSTON EXCELSIOR in the Service of Industry 


Marston’s have unrivalled experience in the fabrication of light alloys—of all shapes and sizes. 


Their products have earned a reputation for efficiency and reliability that is world-wide. 


* Light Alloy Fabrication 

* Specialised Engineering Assemblies 
* Laminated Plastic Components 

* Flexible Tanks 

* Radiators and Heat Exchangers 


Tubular Heat Exchangers 


5 pass cross contra type made tor 
1.C.1. Billingham Division. Length 
18! feet. Diam. 234” inside. Header 
tanks are Aluminium Castings 
of BS.1490 LM.20 alloy 2” thick. 
Body rolled and welded from 
PA.19 Aluminium Alloy Plate {” 
thick. 


Tested to 7§ p.s.i. 


MARSTON EXGELSIOR LIMITED 


(A subsidiary company of Imperial Chemical Industries Ltd.) 
FORDHOUSES, WOLVERHAMPTON. Tel: Fordhouses 3361 


Oil Tank 


of gallon capacity, made 
for I.C.I. Paints Division. 
Length 14 feet. Diameter 7 feet. 
Fabricated from Kynal P3. 
Aluminium alloy {” thick with 
dished ends and a rolled plate 
body. Argon Arc welded 
throughout. Tested to withstand 
an 18 foot head of water. 
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**Newallastic’’ bolts and studs have qualities which are absolutely unique. 
They have been tested by every known device, and have been proved to 
be stronger and more resistant to fatigue than bolts or studs made by 


the usual method. 


/ 
/ 
POSSiILPARR GLASGOW: N 


Atmospheric and Vacuum Distillation Units 


Combined Distillation, Cracking, Reforming and 
Vapour Phase Treating Units 


Pressure Distillate Re-run Units 

Gasoline Recovery and Stabilization Units 
Fractionating Columns and Tube Stills 
Wax Refining, Sweating and Moulding 


A. F. CRAIG & COMPANY LIMITED 


CALEDONIA ENGINEERING WORKS + PAISLEY + SCOTLAND 
LONDON OFFICE: 727 SALISBURY HOUSE - LONDON WALL - £.C.2 + PHONE NATIONAL 3964 
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